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Abstract 



Cosmological perturbation theory is crucial for our understanding of the universe. The 
hnear theory has been well understood for some time, however developing and applying 
the theory beyond linear order is currently at the forefront of research in theoretical 
cosmology. 

This thesis studies the applications of perturbation theory to cosmology and, specif- 
ically, to the early universe. Starting with some background material introducing the 
well-tested 'standard model' of cosmology, we move on to develop the formalism for per- 
turbation theory up to second order giving evolution equations for all types of scalar, 
vector and tensor perturbations, both in gauge dependent and gauge invariant form. 
We then move on to the main result of the thesis, showing that, at second order in 
perturbation theory, vorticity is sourced by a coupling term quadratic in energy density 
and entropy perturbations. This source term implies a qualitative difference to linear 
order. Thus, while at linear order vorticity decays with the expansion of the universe, 
the same is not true at higher orders. This will have important implications on future 
measurements of the polarisation of the Cosmic Microwave Background, and could give 
rise to the generation of a primordial seed magnetic field. Having derived this qual- 
itative result, we then estimate the scale dependence and magnitude of the vorticity 
power spectrum, finding, for simple power law inputs a small, blue spectrum. 

The final part of this thesis concerns higher order perturbation theory, deriving, for 
the first time, the metric tensor, gauge transformation rules and governing equations for 
fully general third order perturbations. We close with a discussion of natural extensions 
to this work and other possible ideas for off-shooting projects in this continually growing 
field. 
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1. Introduction 



Over the last few decades cosmology has moved from a mainly theoretical discipline to 
one in which data is of increasing importance. This is promising, since it means that 
we are no longer confined to the theorists' playground, but instead have observational 
data with which to constrain our models. 

At present the main observable that we have with which to test our theories is the 
Cosmic Microwave Background (CMB) radiation. This is radiation that was produced 
when the universe was around 380,000 years old and had cooled enough to allow elec- 
trons and protons to combine to produce Hydrogen atoms. Its first detection in 1964 
by Penzias and Wilson was hailed as a firm success of the hot Big Bang cosmologi- 
cal model, and experiments have been performed in the years after in order to obtain 
more details of this radiation. The Cosmic Background Explorer (CORE) [23] and 
the Wilkinson Microwave Anisotropy Probe (WMAP) [ ] satellites have since probed 
the anisotropics of the CMB, finding that it is extremely isotropic (up to one part in 
100,000), and the Planck ['] satellite is currently taking data to further increase our 
wealth of data on the CMB. This observation of small anisotropics is very much in 
agreement with the theory, which states that quantum perturbations in the field driv- 
ing inflation produce small primordial density fluctuations which are then amplified 
through gravitational instability to form the structure that exists in the universe today. 

In order to study the theoretical framework of the standard cosmological model (see, 
e.g. Ref. [107], for a particularly lucid review), one uses cosmological perturbation 
theory, which is the main topic of this thesis. The basic idea is quite simple: we model 
the universe as a homogeneous 'background' which has inhomogeneous perturbations 
on top. The perturbations can then be split up order-by-order, with each order being 
smaller than the one before. 

Early studies of the linear order theory were mainly done by the following authors. 
Lifshitz pioneered the early work on perturbations in Ref. [8'.)], which was later extended 
in Ref. [90], with Bonnor considering density perturbations in Ref. [25]. This early work 
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by Lifshitz was conducted in the synchronous gauge, which has since been shown to 
exhibit gauge artefacts if one is not careful [ : ]. Though the authors used intricate 
geometrical arguments in order to remove these unphysical modes, the fact that the cal- 
culations cannot be done easily was far from ideal. Thus, Hawking [63] and Olson [124] 
took a different approach to this problem and attempted the first fully covariant study 
of cosmological perturbations, focussing not on perturbations of the metric tensor, but 
instead on perturbations of the curvature tensor. However, the most pioneering work 
in modern cosmological perturbation theory was completed by Bardeen. In Ref. [16], 
Bardeen presented a systematic method for removing the gauge artefacts by construct- 
ing gauge invariant variables. His work focussed on the two metric potentials \E' and 
$, which correspond to the two gauge invariant scalar metric perturbations in the lon- 
gitudinal gauge. This work was then followed by the two review articles by Kodama 
and Sasaki [ ] and Mukhanov, Feldman and Brandenberger [12!)]. These three articles 
together arguably form the basis of linear metric cosmological perturbation theory. 

An alternative approach to metric perturbation theory is often called the covariant 
approach. The approach defines gauge invariant variables using the Stewart- Walker 
lemma (a perturbation which vanishes in the background is gauge invariant [] 15]) and 
was mostly pioneered by Ellis and collaborators [19-51]. An interesting paper made 
a first step towards relating the covariant approach to the metric approach which was 
written by Bruni et al. [ ]. 

Cosmological perturbations can be decomposed into scalar, vector and tensor pertur- 
bations, as we will show in Chapter 2 and, at linear order, the three types of perturba- 
tions decouple from one another. The scalar modes are related to density perturbations, 
vector modes are vortical or rotational perturbations, and the tensor modes are related 
to gravitational waves. The CMB radiation is polarised, and the different types of 
perturbations induce different polarisations. Scalar modes produce only E-mode (or 
curl-free) polarisation and tensor modes produce only B-mode (divergence-free) polari- 
sation. Vectors produce both, but are usually deemed negligible, since any produced in 
the early universe will be inflated away, or will decay with the expansion of the universe 
[75]. 

However, once we go beyond linear order, different types of perturbation no longer 
decouple and so, for example, vector and tensor perturbations are sourced by scalar 
modes. This mathematical difference between linear and higher orders therefore plays 
an important role in the theory and can result in qualitatively different physics beyond 
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linear order which will, in turn, generate different observational signatures. This is, 
really, the main reason for extending perturbation theory beyond linear order, which 
has been studied by many authors in the last few years [■>, 18, 30, 41, 105, 110, 119, 121- 
123, 148] (see Ref. [11-] for a recent review and a comprehensive list of references on 
second order cosmological perturbations). Perturbation theory beyond linear order is 
the central theme of this thesis. 

This thesis is organised as follows. In the remainder of this introductory Chapter we 
present the standard model of cosmology in a more detailed sense, briefly introducing 
inflationary cosmology, and restating our notation that will be used for the remainder 
of this thesis at the end of the Chapter. In Chapter 2 we introduce the theory of non- 
linear cosmological perturbations up to second order. We present the perturbed metric 
tensor, and energy momentum tensor for a perfect fluid (i.e. a fluid with no anisotropic 
stress) and a scalar field. We consider next the transformation behaviour of the dif- 
ferent perturbations under a gauge transformation, using these to choose gauges and 
define gauge invariant variables. Next, we present a discussion of the thermodynamics 
of a perfect fluid, discussing the pressure and energy density perturbations and the 
definition of the non-adiabatic pressure perturbation, which will play a central role in 
the following work. Finally, to close Chapter 2, we briefly consider how non-adiabatic 
pressure perturbations can arise naturally in multiple fluid or multi-field inflationary 
models. 

In Chapter 3 we continue presenting the foundations of cosmological perturbation 
theory and present the dynamic and constraint equations up to second order. Starting 
with the linear order theory we present the governing equations for scalar, vector and 
tensor perturbations of a perfect fluid without fixing a gauge. We then fix a gauge, 
giving the equations in terms of gauge invariant variables for three different gauges: 
the uniform density, uniform curvature and longitudinal gauges, solving the equations 
for the latter two. We then present the Klein-Gordon equation for a scalar field, and 
highlight the important difference between the adiabatic sound speed and the speed 
with which perturbations travel for a scalar field system. Finally, we investigate the 
perturbations of a system containing both dark energy and dark matter. Having laid the 
foundations with the linear theory, we move on to the second order theory, presenting 
the governing equations for a perfect fluid coming form energy momentum conservation 
and the Einstein equations in gauge dependent form. We then present all equations in 
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the uniform curvature gauge which we will use in Chapter 4, including now only the 
canonical Klein-Gordon equation. Finally, in order to connect with the literature, we 
give the equations for scalars in the Poisson gauge. 

Having now developed the tools for second order perturbation theory, in Chapter 4 
we use the qualitative differences between the linear theory and higher order theory to 
show that, at second order in perturbation theory, vorticity is sourced by a coupling 
between first order energy density and entropy perturbations. This is analogous to 
the case of classical fluid mechanics and generalises Crocco's theorem to an expanding 
background. To show this, we start by defining the vorticity tensor in general relativity, 
and then calculate the vorticity tensor at linear and second order using the fluid four 
velocity and the metric tensor defined in Chapter 2. We then compute the evolution 
equation for the vorticity, making use of the governing equations in Chapter 3. At linear 
order vorticity is not sourced, however there exists a non-zero source term at second 
order when allowing for fluid with a general equation of state that depends upon both 
the energy and the entropy. Having derived this qualitative result, we then give a first 
quantitative solution, estimating the magnitude and scale dependence of the induced 
vorticity using simple input power spectra: the energy density derived in Chapter 3, 
and using a simple ansatz for the non-adiabatic pressure perturbation. 

In Chapter 5 we extend the formalism from the second order theory to third order, 
presenting the gauge transformation rules and constructing gauge invariant variables. 
Then, considering perfect fluids and including all types of perturbation, we present the 
energy and momentum conservation equations and give components of the Einstein 
tensor up to third order. We also give the Klein-Gordon equation for a scalar field 
minimally coupled to gravity. Finally, to close, we conclude in Chapter 6 and present 
possible directions in which one can extend the work presented in this thesis. 

1.1. Standard Cosmology 

We now introduce some elements of standard cosmology in a more quantitative sense 
and, in doing so, define our notation. The basic starting point in cosmology is the 
cosmological principle which states that, on large enough scales, the universe is both 
isotropic and homogeneous. In general relativity, geometry is encoded in the metric ten- 
sor Qn^, or the line element ds"^ = g^pdx^dx^ . The general line element for an isotropic 
and homogeneous spacetime, and thus one which obeys the cosmological principle, takes 
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the form^ [87] 

ds'^ = -dt^ + a^{t) 



+ {de"^ + sin^ Odx' 



in spherical coordinates (t; r, 9, x)i where t denotes the coordinate time, a{t) is a function 
that depends only on time, and K denotes the global curvature of the spatial slices, 
where K = 1 denotes a positively curved, or closed, universe, = a flat universe, 
and i^' = — 1 a negatively curved, or open, universe. This is the Friedmann- Robertson- 
Walker (FRW) metric. An alternative way of representing the FRW metric is 

ds^ = -df + a^{t)%jdx'dx^ , (1.2) 

where jij is the metric tensor on spatial hypersurfaces. 

Since current observations are consistent with a flat, K = universe, which is also in 
agreement with inflation, we adopt this choice henceforth and so write the line element 

as 

ds"^ = -dt^ + a^{t)Sijdx'dx^ , (1.3) 

where 6ij denotes the Kronecker delta. 

We note here the importance of the function a{t), called the scale factor, in an ex- 
panding spacetime. We can picture space as a coordinate grid which expands uniformly 
with the increase of time. The comoving distance, x between two points, which is just 
measured by the comoving coordinates, remains constant as the universe expands. The 
physical distance, r, is proportional to the scale factor, 

r = a{t)x, (1.4) 

and so does evolve with time. Thus, an isotropic and homogeneous universe is charac- 
terised not only by its geometry, but also by the evolution of the scale factor [ ' ]. In 
order to quantify the expansion rate, we introduce the Hubble parameter, H(t) defined 
as 

Hit) = ^ = ^, (1.5) 

CL CL 

where an overdot denotes a derivative with respect to coordinate time, which measures 

^Throughout this thesis we use the so-called 'East coast' metric signature ( — h ++) and the positive 
(+ + +) sign convention in the notation of Misner et al. [i IG]. 
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how rapidly the scale factor changes. 

It is often convenient to use, instead of t, the conformal time coordinate rj defined 
through 

dt , , 

-, (1.6) 



;i-7) 



in terms of which the line element (1-3) becomes 



ds'^ = a^irf) — drf + 6ijdx'^dx^ 



In doing this, we have increased the spatial coordinate grid introduced above to a 
coordinate grid over the entire spacetime. We can furthermore define the conformal 
Hubble parameter 

da/dr] a' . , 

= — , (1.8) 

a a 



where we have used a prime to denote a derivative with respect to conformal time. Then, 
the Hubble parameter in coordinate and conformal time are related to one another by 



n = aH 



1.9) 



Having introduced the metric tensor of an FRW universe, we can now go on to discuss 
the dynamical equations. In general relativity the curvature of a given spacetime is 
encoded in the Riemann curvature tensor, defined as [33] 

po po -pa I -pa -pX -pa -pX f-i -t ri\ 

-ri- piPy — p.u,p fil3,y ~^ Xp'- fj.y Xu'- fil3 ^ \^-^^) 

where F^;^ are the Christoffel connection coefficients, defined in terms of the metric 
tensor and its derivatives as 

= i^g'^^igxp,^ + gx-y,p — 9m^) > (i-H) 

where we have introduced the notation gaii,'^ = d^Qa/B- There are two contractions of 
the Riemann tensor that are particularly useful: the Ricci tensor is given by 



(1.12) 
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and the Ricci scalar, which is the contraction of the Ricci tensor 

R = g'"'R^u. (1.13) 

The Riemann tensor obeys the following identity 

V[Ai?HP'^ = 0, (1.14) 

where the square brackets denote anti-symmetrisation over the relevant indices. Eq. (1-14) 
is often called the Bianchi identity. If we introduce the Einstein tensor G^^, defined as 

then the Bianchi identity implies that the divergence of this tensor vanishes identically: 

V^G^, = 0. (1.16) 

The equation of motion in general relativity is the Einstein equation, 

G^u = R^u — ^QfiuR = 8nGT^i, , (1-17) 

where G is Newton's gravitational constant. T^i, is the energy momentum tensor, which 
describes the energy and momentum of the matter content of the spacetime. We take 
the perfect fluid energy momentum tensor which has the following form 

= (po + Po)Mfo)M(o). + PoS^u , (1.18) 

where p and P are the energy density and pressure, and the fluid four velocity, 
satisfying the constraint u^Uf^ = —1. Note that the subscript '0' denotes the value 
of the quantity in the homogeneous and isotropic background: the importance of this 
notation will become apparent in Chapter 2. In addition to the Einstein equation, there 
also exist a set of evolution equations for the matter variables. These are obtained 
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through the covariant conservation of the energy momentum tensor,^ 

V^T^ = 0. (1.19) 

This equation is obtained from the conservation of the Einstein tensor, Eq. (1.16), and 
by using the Einstein equation, Eq, (1.17). 

Now, let us consider the FRW spacetime. Using the fluid four velocity 

<o) = -(l'0)' M(o)^ = -a(l,0). (1.20) 
a 

and the Christoffel symbols for the (flat) FRW spacetime 

Tl, = n, ro=H5.,, r}o = w,, (1.21) 

r|J. = o, r^o = o, r;., = o, (1.22) 

the Einstein equations (1-17) are then 

SttG 

IT 

A-nG 



n^ = —p,a\ (1.23) 



•H' = —{po + 3Po)a\ (1.24) 

where the first equation comes from the 0-0 component and the second from the trace 
of the spatial Einstein equation. Eqs. (1.23) and (1.24) are called the Friedmann and 
acceleration equations, respectively. The acceleration equation can be rewritten, by 
introducing the constant equation of state Pq = wpo, as 

n' = -^il + 3w)poa\ (1.25) 

where w is the equation of state parameter. Energy conservation gives the continuity 
equation 

p', = -3'H{po + Po), (1.26) 

^It should be noted that this is not the only way one can obtain evolution equations. Instead, varying 
the action with respect to the matter fields will result in the same evolution equations except in 
the case where the system cannot be described by an action (e.g., dissipative fluids [ ]). 
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which can also be rewritten as 

p'^ = -m{l + w)po. (1.27) 
This can then be integrated to give 

Po = Po (^) , (1-28) 

where an overbar denotes the value of a quantity today. The scale factor and Hubble 
parameter can then be shown to evolve as 

We now highlight the evolution of the parameters in the different eras of the universe 
in the hot Big Bang model. The first is the radiation era, where the universe is filled 
with a fluid of particles moving at (or close to) the speed of light, such as photons.^ 
The equation of state parameter for such a fluid is w = 1/3, which gives 

Po oc , and a oc t] . (1.30) 

The next epoch is the matter domination era, where the universe is filled with collision- 
less, non-relativistic particles which better models a universe filled with galaxies. This 
matter is called dust, and is well modelled by a pressureless fluid with equation of state 
parameter w = 0. In this era the energy density and scalar factor evolve, respectively, 

as 

Po (X , and accrf. (1-31) 



1.2. Inflationary Cosmology 

Any summary of modern cosmology would not be complete without a short discussion 
on the inflationary paradigm^, which is a crucial part of the standard cosmological 

^We should note that, genericaUy, one would expect the initial radiation era to be violently disordered. 

However, the existence of the inflationary era guarantees that the radiation era is smooth. We will 

introduce inflationary cosmology in the next section. 
^See, e.g., Ref. [^7] for a more complete treatment of the classical Big Bang problems 
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model. Historically, it was introduced in an attempt to solve some outstanding problems 
in the big bang model. These are: 

• the Horizon problem, that CMB radiation coming from areas of the universe that 
were never in causal contact are observed to have the same temperature; 

• the Flatness problem, that in order for the universe to be so close to fiat today, it 
must have started off very close to fiat; 

• the Relic problem, that no topological relics^ are observed, though they are likely 
produced in the early universe. 

In order to solve these problems, inflation was postulated in the early '80s, simulta- 
neously by Guth [{)(.)], Starobinsky [1 Albrecht and Steinhardt [S] and Linde [91, 92]. 
The precise definition of inflation is simple: it is a period during which the universe 
undergoes accelerated expansion, i.e. 

a>0, (1.32) 

where we have switched to coordinate time in the section for clarity. Another definition, 
equivalent to the first is 

^(^]<0. (1.33) 
dt \aH J ^ ' 

This definition is more intuitive, since 1/aH is the Hubble horizon size, and so inflation 
is defined as a period in which the Hubble size is decreasing. This is precisely the 
condition required to solve the flatness problem. Finally, a third equivalent deflnition 
of inflation is 

Po + 3Po<0, (1.34) 

which, since po is always positive from the weak energy condition [. >_:] , implies a negative 
pressure during inflation. 

In order to be entirely critical of inflation, we should state that inflation does not 
completely remove the initial flatness problem. To begin an inflationary phase, there 
must exist a Planck-scale patch of spacetime that is roughly smooth at an early time. It 



^Topological relics, such as magnetic monopoles are generically produced if the symmetry of a Grand 
Unified Theory is restored in the early universe and then broken spontaneously. Such an abundance 
of relics is higher than observation allows. See, e.g., Refs. [64, f -50] for more details. 
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is hoped that the tuning necessary to achieve this is less than that required to suppress 
the spatial curvature at early times in the absence of an inflationary phase, but this is 
not necessarily true; it may be much worse. Similarly, one might argue that inflation 
does not entirely solve the relic problem, depending on the model of inflation, since some 
hybrid models may produce topological defects as a by-product of their operation. Thus, 
it is too simplistic to claim that inflation solves all the problems completely, however 
any further investigation into this is beyond the remit of this thesis. 

The most popular type of matter which can drive inflation is a scalar fleld. A scalar 
fleld has a Lagrangian density of the form 

£ = p(X,y;), (1.35) 

where X is the relativistic kinetic energy, X = ^g^^ ip ^^ip which, for a homogeneous 
fleld, is then X = — |(/3o^, where v^o(^) is the scalar fleld. Inflationary models can be 
classifled according to the form of the Lagrangian, and whether they contain a single 
scalar fleld or multiple flelds. The most simple single fleld inflationary model has a 
Lagrangian 

p{X,<fo)=X-UM, (1.36) 

where U{(fo) is the potential of the fleld. This simple model can take very different 
forms depending on the choice of this potential function. More exotic models have 
been considered more recently which modify the Lagrangian through changing the de- 
pendence on X as well or instead of changing the potential. These models have been 
of interest since they enable us to evade observational bounds placed on the simple, 
canonical, single fleld models. 

For the simple, canonical model, then, we can write down a pressure and energy 
density for the scalar fleld as 

Po = ^¥30^ + Uiifo) , (1.37) 
Po = 1^0^ - Uiipo) . (1.38) 

Then, using the energy conservation equation (1.26) we obtain the Klein-Gordon equa- 
tion for the homogeneous fleld 



00 + ^T-Lifo + f/,^ = , 



(1.39) 
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and the Friedmann equation (1.23) becomes 



3 V ''^"^ 2 

A useful approximation exists, which consists of neglecting the first term of Eq. (1.39) 
and the last term of Eq. (1.40), to give 

H'^fu, (1.41) 
m^o--U,^, (1.42) 

which, for this approximation to be true, demands two parameters defined as 



e{^,)=AnGy-^\ , (1.43) 

^M = 87rG^, (1.44) 

to be small (i.e. much less than 1). This approximation is called the slow-roll approxi- 
mation, and if it holds guarantees that inflation will occur. In fact, inflation ends when 
e becomes 1. We do not discuss specific models of inflation here, instead pointing the 
interested reader to one of the many textbooks available on the topic [87]. 

To close this brief section on inflationary cosmology, we consider some observational 
signatures of the early universe. Since this is somewhat beyond the main aim of this 
thesis, we simply present some results, and refrain from any derivations. To begin, we 
use the two point correlator of the linear scalar field perturbation, {S(pi6(fi), defined 
later on, which one can then relate to the spectrum of the curvature perturbation that 
will source the anisotropies in the CMB, e.g. the comoving curvature perturbation 
(JZiR-i), which is defined in Eq. (2.143). This can then be evolved forwards using a 
Boltzmann code [86, 138, 154], to give predictions for the anisotropies in the CMB. 

In the notation of the WMAP team [77], the primordial spectrum is then taken to 
be a power law with amplitude A7^(A;o)^ and spectral index 
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where A7^(fco)2 = 2.38 xlQ-^ and ns = 0.969, at the pivot scale of fco = 0.002Mpc~^ The 
observations are compatible with the mostly adiabatic and Gaussian inflationary initial 
condition, but are incompatible with other a priori equally motivated suggestions, such 
as density perturbations induced by topological defects. Constraints on inflationary 
model building are discussed in detail in, e.g., Refs. [■), t*]. Similarly, we can write the 
spectrum for tensor perturbations as 

and we can then define the tensor-scalar ratio, r, as [77] 



(1.47) 



Inflationary models can then be tested according to their predictions for these ob- 
servables. The scalar spectral index is well constrained by WMAP observations, and 
in fact a generic, successful feature of inflation is its ability to generate a near scale 
invariant spectrum [87]. However, the tensor-scalar ratio is less constrained, and dif- 
ferent inflationary models make varying predictions for r (e.g. Refs. [1, 5, 99, 142] and 
references therein). It is hoped that we will be able to narrow the error bars on the 
tensor-scalar ratio by future observations of the polarisation of the CMB from the sky 
(Refs. [1, 21]) and from the ground (e.g. Ref. [n.i]) which will in turn enable us to rule 
out some models of inflation. 



1.3. Notation 

To close the introduction, we briefly state some of the notational conventions that will 
be used in this thesis. 

• We use the mostly positive metric signature, (—+++), and the (+++) convention 
in the notation of [116]. 

• Coordinate time is denoted with t and an overdot denotes a derivative with respect 
to coordinate time; the Hubble parameter is H = a/a. 

• Conformal time is denoted by rj and a prime denotes a derivative with respect to 
conformal time; the conformal Hubble parameter is "H = a' /a. 



1.3: Notation 



23 



• Greek indices {/i, i', ■ ■ ■} cover the entire spacetime and take the range {0 . . . 3}. 

• Latin indices cover the spatial shce and take the range {1 ... 3}. 

• The index (as in nP) denotes conformal time, and an index t (as in ti*) denotes 
coordinate time. 

• The order in the perturbative expansion is denoted with a subscript, 0i, or with 
a subscript enclosed in parentheses if the meaning could be ambiguous, as in u'^^y 

• A comma denotes a partial derivative so, e.g., 

dX dX 
^''^^^'^^^9^' or'WhenX^X(r), X,y^—. 

• A semicolon denotes a covariant derivative with respect to the full spacetime 
metric, g^i,, i.e., 

• The Lie derivative along a vector field is denoted and takes the following 
forms for a scalar 99, a vector, f^, and a tensor, t^^,: 

£^ip = eV,^, (1.48) 
£iv^ = v^,^C + VcC,^., (1-49) 

£itij,,y = t^i/,AC'^ + + tXu^^,H ■ (1.50) 



2. Cosmological Perturbations 



As discussed in the introduction, cosmological perturbation theory is an extremely 
useful, and successful, tool to study the universe in which we live. In this Chapter 
we introduce cosmological perturbation theory in the metric approach (i.e. where we 
consider perturbations to the metric a la Bardeen [l(i]) in a more formal and quantitative 
manner. We introduce the line element for the most general perturbations to FRW, 
and then define the perturbed energy momentum tensor for both a perfect fluid and for 
a scalar field (with both a canonical and non-canonical action). We then derive gauge 
transformations for the different types of perturbation (scalar, vector and tensor), and 
use these to define choices of gauge and gauge invariant variables. We focus on gauge 
choice at linear order and give an illustrative example of how to choose a gauge at second 
order. We close this Chapter by considering the thermodynamics of a perfect fluid, 
defining what we mean by non-adiabatic, or entropic, perturbations, and briefly discuss 
how such perturbations can naturally arise in multi-field or multiple fluid systems. 



We consider the most general perturbations to the flat FRW metric, which gives a line 
element of the form 



We choose a flat background metric because it agrees with observations and is math- 
ematically easier to work with, but should note that all the techniques used in this 
Chapter are valid for a background spacetime with non-zero curvature. 

The perturbations of the spatial components of the metric can be further decomposed 



2.1. Metric Tensor 



ds^ = a^{ri) - (1 + 20)rfr/2 + 2Bidx'dri + {Sij + 2Cij)dx'dx^ 
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as 

Bi = B^i — Si , (2.2) 
Cij = -ip6ij + E^ij + F(^ij) + . (2.3) 

The perturbations are classified as scalar, vector and tensor perturbations according 
to their transformation behaviour on spatial three hypersurfaces [112, 120]. The scalar 
metric perturbations are 0, the lapse function, ip, the curvature perturbation and E and 
B, which make up the scalar shear. Si and Fi are divergence free vector perturbations, 
and hij is a transverse, traceless tensor perturbation. The perturbations therefore obey 
the following relations 

diS' = Q, (2.4) 
d^F' = , (2.5) 
d,h'^ = = h\. (2.6) 

Since the perturbations are inhomogeneous, they depend upon both space and time, 
e.g. 

= 0(x^) = 0(r/,a;^). (2.7) 

The scalar perturbations each contribute one degree of freedom to the perturbed metric 
tensor, each divergence free vector perturbation has two degrees of freedom, as does the 
transverse, traceless tensor perturbation, so we see that, in total, there are 10 degrees of 
freedom - the same as the number of independent components of the perturbed metric 
tensor. Each perturbation can then be expanded in a series: For example the lapse 
function is split as 



UtT^^ (2-8) 



= f^'i + ^tV2+^«''#'3 + -", (2.9) 

where the subscript denotes the order of the perturbation and e is a fiducial expansion 
parameter. We will often omit e when not required for brevity and write 

= 01 + ^02+^03 + -- - , (2.10) 
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In order to define this expansion uniquely, we choose the first order quantity, 0i, to have 
Gaussian statistics. The series is then truncated at the required order. Performing this 
spUt then gives us the covariant components of the metric tensor up to second order 

(700 = -a' (1 + 201 + (62) , (2.11) 
goi = i2Bu + B2i) , (2.12) 
Qij = {Sij + 2Cuj + C2ij) . (2.13) 

The contravariant components of the metric tensor are obtained by imposing the con- 
straint 

g,.g''' = 6,\ (2.14) 
to the appropriate order. To second order this gives 

= (1 - 201 - 02 + 401^ _ B,,B^^ , (2.15) 
gOi = ^(^Bl + - 20ii?i - 2i?ifeCf ) , (2.16) 
= 1 (^S'^ - 2Cf - Ci' + 4Cf Cife^' - BiB^j . (2.17) 

2.2. Energy Momentum Tensor 

The matter content of the universe is described by the energy momentum tensor. Since 
General Relativity links the geometry of spacetime to its matter content, perturbations 
in the metric tensor invoke perturbations in the energy momentum tensor.^ In this 
section, we outline the perturbed energy momentum tensor for a perfect fluid and a 
scalar field, to second order in perturbation theory. 



There are, in fact, works where this is not the case and, for example, only the energy momentum 
tensor is perturbed. However, in order for the work to be consistent, one should perturb both the 
geometry and the matter content of the spacetime. 
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2.2.1. Perfect Fluid 

The energy momentum tensor for a perfect fluid, i.e. in the absence of anisotropic 
stress, is as presented in the previous section, 

= (p + P)u^u, + P5\ , (2.18) 

where p is the energy density, P is the pressure and is the four velocity of the fluid. 
Note that, for the purposes of this thesis, we define a 'perfect fluid' to be a fluid with 
a diagonal energy momentum tensor as in, e.g., Ref. [17]. 

The energy density and the pressure are expanded up to second order in perturbation 
theory in the standard way 

P = Po + 5pi + ^5p2, (2.19) 
P = Po + 5Pi + ^5P2 . (2.20) 

The fluid four velocity, which is defined as 

= (2.21) 

where r is an afline parameter, here the proper time, and is subject to the constraint 

u% = -l. (2.22) 

To second order in perturbation theory the fluid four velocity has the contravariant 
components 

u'=^- {v,' + ^^2^^ , (2.23) 
M° = ^ (^1 - 01 - ^02 + \<Pl + + viuBi^^ , (2.24) 
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and the covariant components, 

Ui = a (^u + Bu + ]^{v2i + B2i) - <PiBu + 2CukVi'^ , (2.25) 

Mo = -a (^1 + 01 + ^02 - ^0? + ^t^i'^t^ifc^ , (2.26) 

where f * is the spatial three velocity of the fluid. Then, by substituting the components 
of the four velocity along with the expansions of the energy density and pressure into 
Eq. (1.18), we obtain the components of the energy momentum tensor, up to second 
order 

T% = -(Po + Po)(^i' + 5i'=)t;i, - (po + hi + \h2) , (2.27) 
T°, = (po + Po) (^u + Bu + ]^{v2i + B2i) - <pi{vu + 2Bu) + 2CukVi' 

+ {5pi + 5Pi){vu + Bu), (2.28) 
T'j = {P, + 5Pi + ^5^2^ + (po + Po)^i'(^i, + 5i,) . (2.29) 

2.2.2. Scalar Field 

As shown in Section 1.2, scalar fields play an important role in modern cosmology 
through the theory of inflation. The energy momentum tensor for a scalar field min- 
imally coupled to gravity is defined through the Lagrangian density. For a canonical 
scalar field, the Lagrangian takes the form 



C = --g^''v,^V,.-U{^), (2.30) 



where U{ip) is the potential energy of the scalar field. The variational energy momentum 
tensor (or Hilbert stress-energy tensor) is then defined as [87] 

dC 

T,^^-2— + g,X, (2.31) 
and, for a scalar field if, we obtain 

= g''^V,xV,, - 5^u {u{^) + ]^g^^^^^^^p\ . (2.32) 
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Expanding the scalar field in the usual way, and using the definition of the metric tensor 
given in Section 2.1, gives the components of the energy momentum tensor for a scalar 
field up to second order in perturbation theory (see, e.g., Refs. [3, 105, 123]) 



2a2 L 



1 - 201 + 2^ + 40i2 - B^kB^, + 20i Vi - f ^ 



+ 



V2 



TV 



51 + ^ - 20151 - 2B,,C- + 2^ + '^'^'^^^ 
-2^^ + ife + ^^^-' 



(2.33) 



2 5 



(2.34) 



L ' (y^Q ' 2 ' 



(2.35) 



2a2 



1 - 201 + 2M + 401^ - 5i.5f + 4^ - 

^0 "Po "Po 



^0 



V2 



(2.36) 



We can write down a pressure and energy density for the scalar field by comparing 
the components of the scalar field energy momentum tensor to that of a perfect fiuid 
given in Eqs. (2.27)-(2.29). In the background this gives, as shown in Section 1.2, 



1 

2^2 



Pn 



2a2 



(2.37) 



and, to linear order, 



(2.38) 
(2.39) 



We can also express the fiuid velocity in terms of the field by comparing Eq. (2.28) to 
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Eq. (2.35) to give 

= (2.40) 

where Vi = Bi + Vi. 

There has been a lot of recent interest in scalar fields with non-canonical actions: 
early work in the realm of the early universe such as k-inflation [12, 56] and more 
current work including the string theory motivated Dirac-Born-Infeld (DBI) infiation 
[9, 66, 88, 135, 140] . Scalar fields with non-canonical actions have also recently been 
considered as dark energy candidates (see, e.g. Ref. [149]). These models all modify 
the Lagrangian (2.30) to the more general function^ 

C = p{X,^), (2.41) 

where X = —^g^'^ip^^ip ^. The Lagrangian for DBI infiation is then given by 

p(X, ^) = -T(^) v/l-2T-i(^)X + T(^) - V{^) , (2.42) 

where T{ip) and V{ip) are functions that further specify the model. 

The energy momentum tensor for the general Lagrangian (2.41) is then obtained from 
Eq. (2.31) as 

T% = p,x9^^^,x^,, + . (2.43) 

This has the components, up to linear order in perturbation theory, and in coordinate 
time, 

T\ = {po - 2p,xXo) - [{p,x + 2Xop,xx) SXi + 2X^p^x^5^i - , (2.44) 

T^ = -p,x¥3o5<^i,., (2.45) 
T', = (po + P,xSX, + P,^5ip,) 6', , (2.46) 

where we have defined SXi = ipoSip^ — 0i</3o^, and Xq = |(/3o^- As in the canonical 
case, we can define a pressure and energy density for the scalar field. The background 
quantities are 

Po = 2p,xXo-p, Po = P, (2.47) 

^One could, of course, cook up even more general Lagrangians. However, typically if higher-order 
derivatives of the scalar field are present then the model will exhibit Ostrogradski instabilities (see 
Ref. [47] for an outline of Ostrogradski's original argument). 
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and the linear order quantities are 



5Pi =p^x5Xi+p^^5(f^ . 



(2.48) 
(2.49) 



2.3. Gauge Transformations at First and Second 



Gauge transformations play an important role in cosmological perturbation theory. 
General relativity is a theory of differential manifolds with no preferred coordinate 
charts, and is therefore required to be covariant under coordinate transformations. 
However, when we come to consider perturbations in general relativity we must, for 
consistency, consider perturbations of the spacetime itself. That is, in the language 
of differential geometry, we consider a one parameter family of four-manifolds M^, 
embedded in a five-manifold N [108, 144]. Each represents a spacetime, with the 
base spacetime, or unperturbed e = manifold, Mq. The problem comes in perturbation 
theory when attempting to compare two objects which 'live' in different spaces. In order 
to deal with this, we introduce a point identification map p^ : Mq — )■ which relates 
points in the perturbed manifold with those in the background. This correspondence 
introduces a new vector field, X on A^, and points which lie on the same integral curve 
7 of X are regarded as being the same physical point. 

However, the choice of this point identification map and, therefore, the vector field 
X is not unique. The choice of the correspondence between the points on the Mq and 
those on or, equivalently, the choice of the vector field X is called a choice of gauge, 
and X is then the gauge generator. A gauge transformation then tells us how we move 
from one choice of gauge to another. 

There are two approaches to gauge transformations. First, consider a point p on Mq. 
Two generating vectors X and Y then define a correspondence between this point p 
and two different points s and t on M^. Clearly, then, these choices induce a coordinate 
change (gauge transformation) on M^. This is known as the passive view. Alterna- 
tively, consider a point p on M^. We then find a point s on Mq which maps to p under 
the gauge choice X and a point t, also on Mq that maps to p under the choice Y. In 
this case, a gauge transformation is induced on the background manifold, Mq. This is 



Order 
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known as the active view. We can think of the active approach as the one in which the 
transformation of the perturbed quantities is evaluated at the same coordinate point, 
and the passive approach where the transformation of the perturbed quantities is taken 
at the same physical point. 

In this section we go on to discuss the active and passive approach briefly then, 
adopting the active approach, we derive the gauge transformation rules for scalars, 
vectors and tensors up to second order in the perturbations. 



2.3.1. Active Approach 

In the active approach to gauge transformations, the exponential map is the starting 
point [108, 119]. Once the generating vector of gauge transformation, has been 
specified, we can immediately write down how a general tensor T transforms. The 
exponential map is 

T = e^«T, (2.50) 

where £^ denotes the Lie derivative with respect to the generating vector which, up 
to second order in perturbation theory, is 

e ^ + \eH'2 +■■■■ (2.51) 

The exponential map is then 

e-K^{£i) = 1 + e£^, + ^e^£l + ^e^£^, + ■■■ (2.52) 

up to second order in perturbation theory. Splitting T order by order, we find that the 
tensorial quantities transform at zeroth, first and second order, respectively, as [30, 119] 

To = To, (2.53) 
eSTi = eSTi + e£^,To , (2.54) 
e^6T2 = (^Ta + ^^.Tq + ^^Tq + 2£^,6Ti^ . (2.55) 



By noting that the Lie derivative acting on a scalar is just the directional derivative, 
£^ = ^^{d/dx'^), the exponential map can be applied to the coordinates to obtain 
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the following relationship between coordinates at two points, p and 



x^{q) = e 



ip)- 



(2.56) 



Expanding this to second order gives 



(2.57) 



This coordinate relationship is not required to perform calculations in the active ap- 
proach, but will be useful for the discussion in Section 2.3.2 below. 

2.3.2. Passive Approach 

A natural starting point for discussing the passive approach to gauge transformations 
is the coordinate relationship Eq. (2.57) since, in the passive approach, one states the 
relationship between two coordinate systems and then calculates how variables change 
when transforming from one coordinate system to the other. However, since in the 
passive approach quantities are evaluated at the same physical point we need to rewrite 
Eq. (2.57) Choosing p and q such that xf^{q) = x'^{p), Eq. (2.57) enables us to 

write 



a;^(g) = x^{p) 



(2.58) 



Using the first terms of Eq. (2.57), 



x^{q)=x^{p) + eeAp). 



(2.59) 



allows us to Taylor expand as 



= er(g) - <.(g)ei\g) , 



(2.60) 



■^Note that, in some of the literature, a different sign is taken in the exponent in this equation to 
obtain correspondence between the active and passive approach at hnear order. Since we are not 
solely working at linear order in this work, there is no advantage gained by making such a choice. 
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which is vahd up to second order. Substituting Eq. (2.60) into Eq. (2.58) gives the 
relationship between the two coordinate systems at the same point, q 



xt^{q) = x^iq) - eCf (g) 



(2.61) 



Having highhghted the two approaches to gauge transformations we now focus on 
the active approach and give some concrete examples. 



2.3.3. Four Scalars 

We will now look at the transformation of four scalars, choosing the energy density, p, 
which can be expanded as in Eq. (2.19), 

p = Po + 6pi + hp2 + ^Sps + ■■■ , (2.62) 

as an example. 



First Order 

Before studying the transformation behaviour of perturbations at first order, we split the 
generating vector C,i into a scalar temporal part ai and a spatial scalar and divergence 
free vector part, respectively f3i and 71*, as 

er = (ai,/3i; + 7ir (2.63) 
Since the Lie derivative of a scalar p with respect to the vector ^'^ is simply 

£^p = ep,,, (2.64) 
from Eq. (2.54), we then find that the energy density transforms, at linear order, as 

Spi = Spi + p'^ai . (2.65) 

We see that, at first order, the gauge transformation is completely determined by the 
time slicing, ai.^ 

^Note that ai does not generate a foliation of spacetime by spatial hypersurfaces - this is inherited 
from the foliation already present in the background spacetime. Instead, ai labels the time slicing. 
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Second Order 

At second order we split the generating vector ^2 analogous way to first order as 

e2^ = (a2,/32;+72r (2.66) 

Then, using Eq. (2.55), we find that the second order energy density perturbation 
transforms as 

5p2 = ^P2 + Poa2 + ai(Po«i + Poa'i + 25pi) + {26pi + p'^^ai) ^k{(^i^ + li") ■ (2.67) 

Thus, at second order, the gauge transformation is only fully determined once the time 
slicing is specified at first and second order (ai and 0:2) and the spatial threading (or 
spatial gauge perturbation) is specified at first order (/3i and 71*) [112]. 

2.3.4. The Metric Tensor 

We will now focus on the transformation behaviour of the metric tensor. Again, the 
starting point is the Lie derivative, which for a the metric tensor, gf^jy, is given by 

£^g^u = g^.u,xt + g^.xi\ + ^^a.^^,^ ■ (2.68) 

First Order 

At first order, the metric tensor transforms, from Eqs. (2.54) and (2.68) as 

^} = Sg^J + g':U^ + g':k^ + g'S^^ „ . (2.69) 

We can obtain the transformation behaviour of each particular metric function by 
extracting it, in turn, from the above general expression using the method outlined in 
Ref. [112]. From Eq. (2.68) we obtain the following transformation behaviour for Cuj 

2Cuj = 2Cuj + 2Hai5ij + + • (2.70) 

From this we can extract the transformation behaviour of the spatial metric functions. 
Here we do not focus on the details, but instead quote results. We find that the scalar 
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metric perturbations transform as 



01 = 01 + -Htti + a'l , (2.71) 

V^ = ^i-^ai, (2.72) 

Bl = Bi-ai + l3[, (2.73) 

^ = ^1 + /3i , (2.74) 



the vectors as 



5i^ = 5i^-7i , (2.75) 
Fr* = Fi^ + 7iS (2.76) 

and, as is well known, the tensor component, huj, is gauge invariant. Finally, the scalar 
shear, which is defined as 

ai = E[-Bu (2.77) 

transforms as 

ai = ai + ai, (2.78) 
which will be useful later when we come to define gauges and gauge invariant variables. 

Second Order 

At second order we obtain the transformation behaviour of the metric tensor from 
Eqs. (2.55) and (2.68), noting that 

£|T = %(%T). (2.79) 

The metric tensor therefore transforms as 



I (0) fAfO I (0) f a I O „vu; fOfA , v^; ^a^A , „Vu;fA 

ASl ,aSl I ^ y/iA,QSlSl ,1/ I yA i/,aS>l SI I i/AaSl ,^S1 ,1/ 



(0) t-a^-A I ^(0) (-acX , ^(0)£.A 



+ {ii + ei .) + gf^ i^' + ei ,) , (2.80) 
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from which we can extract, as at first order, the transformation behaviour of individual 
metric perturbation functions. It is a httle trickier to obtain the transformation be- 
haviour of iIj2 from this expression. However, we note that the expression in Eq. (2.80) 
gives the transformation of C2ij, namely,^ 



2C. 



2ij 



(2.81) 



where Xij contains terms quadratic in the first order perturbations and is defined below 
in Eq. (2.87), and so we need intermediate methods in order to extract the transforma- 
tion of a particular component. Again, we we do not go into the unnecessary details 
here, but instead refer the interested reader to Ref. [112]. After this calculation we 
obtain 

= - -Ha^ - \x\ + . (2.82) 



We find that the other second order scalar metric perturbations transform as [112] 



+ Ha2 + a2 + ai [ai + bUai + [W + 2%^) ai + 4H0i + 2^;] (2.83) 



+ 2ai' (ai' + 20i) + ^a. (ai' + Ha^ + 20i) ' + 



ar - 25u, 



^2 
B2 



B2-a2 + f5'2 + V-^A-B^ 



k 1 



(2.84) 
(2.85) 



where X-qi and Xij are defined as 



{2HBu + B[^) «i + Bu,kii - 20i«i,,, + , + B^a'^ + 2Ci,kii 



(2.86) 



^In the following and for the rest of this section, we do not split up the spatial part of the gauge 
transformation generating vector into scalar and vector parts for brevity. We denote the spatial 
part of by 
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and 



Xij = 2 



+ 4 



a 



(2.87) 



Furthermore, the vector perturbations transform as 



S2i — S2i — 72/ — + V '^X'B^^ki 1 



Mi 1 



(2.88) 
(2.89) 



and the tensor perturbation which at second order, unhke at first order, is not gauge 
invariant, as 



1 



1. 



,klij 



(2.90) 



2.3.5. Four Vectors 

Finally, we move on to the transformation behaviour of a four vector. The Lie derivative 
of a vector, PV^, is given by 



£iW^ = w^,^c + w^e 



First Order 

A four vector transforms at first order, from Eq. (2.91), as 

which gives, for the specific case of the four-velocity, the transformation rule 



(2.91) 



(2.92) 



Vii = Vii - ai, 



(2.93) 
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where the quantity Vu is defined as 

Vu = vu + Bu. (2.94) 

Then, on sphtting the velocity perturbation into a vector part and the gradient of a 
scalar, as 

vii = + f i,i , (2.95) 
recalling that the metric perturbation Bn can be split up as 

Bu = - Su , (2.96) 

and making use of Eqs. (2.73) and (2.75), we obtain the transformation rules 

vi = vi-(3[, (2.97) 

vfi = v'^' - 71^' . (2.98) 

Second Order 

At second order we find that a four vector, VF2^j, transforms as 

= + l^('o)^«2 + W^(o)o«2,;. + + W^('o)^«i,,ei' , (2.99) 



A 



Focussing again on the fiuid four velocity, we obtain from the ith component of 
Eq. (2.99), the transformation behaviour of V2i: 

= V2^ - a2,i + 4ef Qfc + 2a[ (i[ + Bu - + ai(2V/, - «;,,) 

+ 2(6.,fc + iikMi - ^i") + '2na,{2Bu - Vu + 3^1,) 

+ (2V^H,fc + «i,^fc) + ^W^V,k - ai,fc) + 201 (el - 2ai,,) . (2.100) 

Then, using the transformation behaviour of i?2i, given by 

B^2i = B2i+C2i-a2,i + XBi, (2.101) 
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we find that the second order three- velocity, V2i, transforms as 

+ ^i^u,k + 2ai« - nvu) - 2v'l^u,k + 2vu,k^'l ■ (2-102) 

2.4. Gauge Choices and Gauge Invariant Variables 

As mentioned previously, a central element to Einstein's theory of general relativity is 
the covariance of the theory under coordinate reparametrisation. However, a problem 
arises when undertaking metric cosmological perturbation theory since the process of 
splitting the spacetime into a background and a perturbation is not a covariant process 
(see, e.g., Ref. [48] and Section 2.2 of Ref. [108]) Therefore, in doing so, one introduces 
spurious gauge modes so that variables depend upon the coordinate choice. Observa- 
tional quantities should not depend upon the choice of coordinate used, and therefore 
this so-called gauge problem of perturbation theory seemingly would introduce confu- 
sion and erroneous results. However, as long as one is careful to remove the gauge 
modes, this will not pose us a problem. 

The gauge problem was first 'solved' by Bardeen in a consistent way in Ref. [!()], and 
has been studied in much detail, and extended beyond linear order, in the decades since. 
The solution lies with the introduction of gauge invariant variables, that is, variables 
which no longer change under a gauge transformation. Bardeen constructed two such 
variables for scalar perturbations, which happen to coincide with the lapse function 
and curvature perturbation in the longitudinal gauge (see Section. 2.4.2). However, 
the systematic approach can be extended to other gauges: one simply inspects the 
gauge transformation rules presented in the previous sections, and chooses coordinates 
such that, e.g., two of the scalar metric perturbations are zero. This enables one to 
remove the gauge dependencies ai and Pi, rendering the other scalar perturbations 
gauge invariant. Similarly, this can be extended beyond scalar perturbations, and one 
can inspect the transformation rules for a vector perturbation, setting it to zero, and 
thus removing the dependency on Yi- 

Finally, let us emphasise the 'gauge issue' by considering degrees of freedom of the 
metric. In four dimensions, the metric starts with 16 degrees of freedom: 6 are lost 
because of symmetry, 4 more because of coordinate invariance [gauge choice) of the 
metric and 4 to do with the Hamiltonian constraints, which arise when writing down 
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the field equations - also called a gauge choice. This leaves 2 degrees of freedom for 
the two polarisations of the graviton. But this is so far not related to perturbation 
theory. By perturbing the metric, one introduces further degrees of freedom which are 
not addressed by the above - the option to change the gauge, or map. These are the 
degrees of freedom addressed by the perturbation gauge choice. 

This paragraph does well to highlight that 'gauge' is a well used term in theoretical 
physics, and is often used to mean different (albeit closely related) things. When dis- 
cussing cosmological perturbation theory, and for the rest of this thesis, we reserve the 
phrase 'choice of gauge' to mean a specification of the mapping between the background 
and the perturbed spacetimes. 

In this section we concentrate on the linear theory, present the definitions of various 
gauges commonly used throughout the literature, and define some gauge invariant vari- 
ables. We then give an example of how the theory works at second order, by describing 
the uniform curvature gauge. This section is mainly a review, and more details can be 
found in Ref. [112]. 

2.4.1. Uniform Curvature Gauge 

A possible choice of gauge is one in which the spatial metric is unperturbed. At linear 
order, this amounts to setting Ei = ipi = and Fu = 0. This specifies the gauge 
generating vector, , using Eqs. (2.72) and (2.74), as 



i^i = i^i — 'Hoii 







(2.103) 



^ «lfiat 



n 



(2.104) 



El = El + ^1 = 

=^ /^Iflat = —El , 



(2.105) 
(2.106) 



and 



Fi^ = Fi' + 7i' = 
^ 7iflat^ = -Fi' 



(2.107) 
(2.108) 
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The other scalars in this gauge, which are then gauge invariant are, from Eqs. (2.71) 
and (2.73), 



B^t = Bi-^-Ei'. (2.110) 
The gauge invariant vector metric perturbation is, from Eq. (2.75), 

w = ^i* + i^i;. (2.111) 

Note that the scalar field perturbation on spatially flat hypersurfaces is the gauge 
invariant Sasaki-Mukhanov variable [118, 133], often denoted by Q, 



2.4.2. Longitudinal (Poisson) Gauge 

The longitudinal gauge is the gauge in which the shear, a, vanishes. It is also known 
as the conformal Newtonian or orthogonal zero-shear gauge. Its extension to include 
vector and tensors is called the Poisson gauge. This gauge is commonly used in the 
literature, since the remaining gauge invariant scalars in this gauge are the variables 
introduced by Bardeen [Ki]. 

At linear order the temporal part of the gauge generating vector is specified by the 
choice oi = as, using Eq. (2.78) 

a^ = cxi + ai = (2.113) 
^au = -ai, (2.114) 

where the subscript i denotes the value in the longitudinal gauge. The generating 
vector is fully specified, for scalar perturbations, by making the gauge choice Ei = 
(and hence Bi = 0) as 

e'i = Ei + (3i=0 (2.115) 
^/3h = -Ei. (2.116) 
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The other two scalar metric perturbations in this gauge are then 



0i, = 0i--Hai-a;, (2.117) 
= ^i+Hai, (2.118) 

which, by using the definition of the shear, Eq. (2.77), give 

= 01 - -HiE,' - B{) - (E/ - B{)' , (2.119) 
iru = iJi+V-{Ei' -B^). (2.120) 

These are then identified with the two Bardeen potentials, $i and \l'i, respectively (or, 
in Bardeen's notation, $aQ^°^ and -$hQ^°^)- 

In the Poisson gauge, the generalisation of the longitudinal gauge beyond scalar 
perturbations, the spatial vector component of the gauge transformation generating 
vector by demanding that S'l* = gives, using Eq. (2.75), 

S^i = Si' + = Q (2.121) 
=^7i/ = j Si'dv + J'lKx^) , (2.122) 

where J-*!* is an arbitrary constant three-vector. Thus, when having fixed the Poisson 
gauge, there still exists some residual freedom in this choice of constant vector. 

The remaining gauge invariant vector metric perturbation in the Poisson gauge is 
then 

= Fi^ + J Si'dT] + J^i\x^) . (2.123) 
2.4.3. Uniform Density Gauge 

As an alternative to the gauges above, we can define a gauge with respect to the matter 
perturbations. One example is the uniform density gauge which is based upon choosing 
a spacetime foliation such that the density perturbation vanishes. 
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At first order we can fix ai by demanding that 6pi = 0. Using Eq. (2.65) we obtain 



We still have the freedom to choose the spatial scalar part of the gauge transformation 
generating vector, which can be done unambiguously by choosing, e.g., Ei = 0. 

An especially interesting variable is the curvature perturbation in this gauge, (i, 
which is a gauge invariant variable and defined as 



where the sign is chosen to agree with that in Ref. [j i]^. This quantity will come in 
useful in later chapters because it is conserved on large scales for an adiabatic system, 
as will be shown in Section 3.1.1. 

2.4.4. Synchronous Gauge 

The synchronous gauge was popular in early work on cosmological perturbation theory, 
and was introduced by Lifshitz in the groundbreaking work of Ref. [89]. It is char- 
acterised by 01 = = Bii, so that the Qqq and goi components of the metric are left 
unperturbed and any perturbation away from FRW is confined to the spatial part of 
the metric. It can be thought of physically as the gauge in which rj defines proper time 
for all comoving observers. This gauge is also used in many modern Boltzmann codes 
such as CMBFAST [138], and is discussed in detail, and compared to the longitudinal 
gauge in Ref. [102]. 

The synchronous gauge condition fixes the temporal gauge function through 



5pi = 5pi + PqOi = 



(2.124) 




(2.125) 




(2.126) 



01 = 01 + "Htti + a[ = 
=^ a0i + a'aisyn + aa[ 



Isyn 



a0i + (aaisyn)' = 



(2.127) 
(2.128) 




(2.129) 



^See Ref. [ ] for a detailed comparison of different sign conventions and notation used for the 
curvature perturbations in different papers. 
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and the spatial gauge functions as 

/3isyn = j{ai,y,,-Bi)dri + Vix'), (2.130) 
7Wn= / Sldv + S\x'). (2.131) 



isyn 

The function V^i{x'') + £i{x^) affects the labelling of the initial spatial hypersurface. 
However, the function C{x^) affects the scalar perturbations, and so the synchronous 
gauge does not determine the time-slicing unambiguously. It is therefore not possible 
to define gauge invariant variables from the metric in this gauge, since the remaining 
scalars (for example the curvature perturbation), have spurious gauge dependence: 

= ^1 + ^ ( y a<i)^dri - C(x'=)) . (2.132) 

In Lifshitz' original work, the gauge mode was removed using symmetry arguments. 
Nowadays a systematic approach is used to remove this gauge mode, and a further 
gauge condition is taken, setting the perturbation in the three velocity of the dark 
matter fluid to zero. Then, 

v^Id^yn = vi- I3[^y^ = (2.133) 
^ a{vi + Bi) + j a(i)idri - C{x^) = . (2.134) 

However, we then refer to the field equations, Eq. (3.2) which, in the synchronous 
gauge guarantees that for the cold dark matter perturbation (where = = 6Pi), 
Vicdm = «,(fi + -Bi)|cdm is a constant (in conformal time). Thus, we obtain 

C{x') = a{vi + Bi) , (2.135) 

cdm 

which removes the gauge freedom. Note that we can only choose the dark matter fluid 
with which to define the gauge since it is pressureless. The same does not hold true 
for a fluid with non-zero pressure. In order to see this we need to use an equation that 
we will derive in full detail later from the momentum conservation, Eq. (3.7) which, in 
coordinate time, is 

Vi - mclVi + + (/>! = . (2.136) 

PO + -TO 
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In the synchronous gauge, i?i = = 0i, so this becomes 

avi - mclavi + ^^^^^ = . (2.137) 

If we then set the three- velocity to zero, in order to specify the threading, Eq. (2.137) 
becomes 

cl5p 



77=0, (2.138) 



which states that the fluid is pressureless. Therefore, we can only define the synchronous 
gauge as comoving with respect to a pressureless fluid, e.g. cold dark matter, and not 
with respect to a fluid with pressure. 

2.4.5. Comoving Gauge 

Another example of a gauge defined by the gauge transformation of a matter variable is 
the comoving gauge. This is defined by choosing the gauge such that the three-velocity 
of the fluid vanishes, Vu = 0. Furthermore, this choice implies that Vu = 0. Then, 

Vi = Vi-ai = (2.139) 

^ aicom = Vi = vi + Bi, (2.140) 

and 

vi = vi-l3[ = (2.141) 

/3icom = j vidr] + J^ix") , (2.142) 

where J^{x^) denotes a residual gauge freedom. However, note that scalar perturbations 
are all independent of J^{x^). One of the remaining scalars is the curvature peturbation 
on comoving hypersurfaces, which is quite popular in the literature, and often denoted 
7^: 

n = = i^i - nvi . (2.143) 
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2.4.6. Beyond Linear Order 

At second order the procedure is very much the same as at hnear order, and the various 
gauges are defined in an analogous way to linear order. Of course, the expressions 
obtained are much longer, due to the fact that the second order gauge transformations 
contain many more terms than those at first order. Since we do not intend this work to 
be a comprehensive summary of gauge choice and gauge invariant variables (this topic 
has already been covered in full, gory detail in Ref. [11-]), we instead sketch how the 
gauge choice and construction of gauge invariant variables work at second order for one 
choice of gauge: the uniform curvature gauge. 

As at linear order, detailed above in Section 2.4.1, we determine the components of 
the gauge transformation generating vector through the conditions ip2 = E2 = and 
= 0. The first gives, from Eq. (2.82), 



flat,«i 



Ml" 



flatfc i 5 



(2.144) 



where Aflatij is defined as Eq. (2.87) with the first order generators given above in 
Eqs. (2.104), (2.106) and (2.108), and the second condition gives 



f3: 



2flat 



2 -yij 



-E2 - -VVX, 



(2.145) 



Finally, imposing the condition F2i = 0, gives 



72flati 



^2i 



flatifc 



(2.146) 



As an example of a gauge invariant variable at second order, we present the lapse 
function in the uniform curvature gauge, 02flat- Using Eq. (2.83) along with Eqs. (2.104) 
and (2.144), we obtain 




4 . . 1 



Ai + 2H0 



/■A; / 
Slflat ; 



(2.147) 
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where we have defined 



A. 



in) 



n 



l^(n) 



(2.148) 



2.5. Thermodynamics of a Perfect Fluid 

Considerable physical insight can be gained by studying the thermodynamic properties 
of a system [80, 96]. In this section, we study a single perfect fluid system. Such a 
system is fully characterised by three variables, of which only two are independent. Here 
we choose the energy density, p and the entropy, as independent variables, with the 
pressure being given by the equation of state P = P{p, S). The pressure perturbation 
can then be expanded, at linear order in perturbation theory, as 



SPi 



dP 

as 



dP 

p op 



5pi 



(2.149) 



This can be cast in the more familiar form 



8Pi = (5Pnadi + cl6pi , 



(2.150) 



by introducing the adiabatic sound speed 



dP 

dp 



(2.151) 



and by defining the non-adiabatic pressure perturbation, which is proportional to the 
perturbation in the entropy, as [152] 



^-Pnadl = 



dP 

'dS 



6Si. 



(2.152) 



Note that (5Pnadi is gauge invariant. This can be shown by considering the gauge 
transformation for the energy density perturbation, Eq. (2.65), along with the analogous 
equation for the pressure perturbation. One can extend Eq. (2.149) to higher order by 
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simply not truncating the expansion at linear order, that is 

dP dP 



1 

+ 2 



52P^^2 d^P , d'^P ^ 2 



+ ... . (2.153) 



The entropy, or non-adiabatic pressure perturbation at second order, for example, is 
then found from Eq. (2.153), as [110] 

<5P2„ad = 5P2 - cl5p2 (2.154) 

op 

2.5.1. Entropy or Non-Adiabatic Perturbations from Inflation 

One way in which a non-adiabatic pressure perturbation can be generated is through the 
relative entropy perturbation between two or more fluids or scalar fields. For example, 
the relative entropy or isocurvature perturbation, at first order, between two fluids 
denoted with subscripts A and B is [^^ ,] (dropping here the subscripts denoting the 
order, for brevity) 

S.s = m{'^-^A. (2.155) 

In a system consisting of multiple fluids, the non-adiabatic pressure perturbation is 
split as [75, 112] 

(5Pnad = 5Pintr + 5Prel, (2.156) 

where the first term is the contribution from the intrinsic entropy perturbation of each 
fluid, and the second term is due to the relative entropy perturbation between each 
fluid. Sab, and is defined as 

^^-1 ^ T^T^PoaPobU - cI)Sab , (2.157) 

where c\ and are the adiabatic sound speed of each fluid. Thus, for a multiple fluid 
system, even when the intrinsic entropy perturbation is zero for each fluid, there is a non- 
vanishing overall non-adiabatic pressure perturbation. This can be extended to the case 
of scalar fields by using standard techniques of treating the fields as fluids. Much recent 
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work has been focussed on the discussion of entropy, or isocurvature perturbations in 
multi-field inflationary models. See, e.g. Refs. [7, 19, 31, 32, 54, 59, 76, 78, 81, 82, 84, 
93, 94, 101, 127, 134] and references therein. 



3. Dynamics and Constraints 



In this Chapter we give the governing equations for perturbations of a FRW universe. 
The background evolution and constraint equations are presented in Chapter 1, so here 
we consider the equations at first and second order in perturbation theory. Starting with 
the hnear order theory we present the governing equations for scalar, vector and tensor 
perturbations for a universe filled with a perfect fluid without fixing a gauge. We then 
make three choices of gauge, the uniform density, uniform curvature and longitudinal 
gauges, and solve the evolution equations for scalar perturbations in the latter two 
case. Next, we present the evolution equation for a scalar field - the Klein-Gordon 
equation - for a field with both a canonical and non-canonical Lagrangian, highlighting 
the importance of the difference between the adiabatic sound speed and the phase 
speed for a scalar field system. We finish our discussion of linear perturbations with 
an investigation into the perturbations of a system with both a dark matter and a 
dynamical dark energy component. 

Having discussed linear order perturbations, we then move on to the second order 
theory. We derive the governing equations for a perfect fluid in a perturbed FRW uni- 
verse from energy momentum conservation and the Einstein equations, without fixing 
gauge. We go on to present the equations for scalar and vector perturbations in the 
uniform curvature gauge, which will come in useful in Chapter 4, and then for scalar 
perturbations only, including now the canonical Klein-Gordon equation at second order. 
Finally, to connect with other parts of the literature, we give the equations for scalars 
in the Poisson gauge. 

3.1. First Order 

In this section we give the evolution and constraint equations at first order in cosmolog- 
ical perturbation theory for a universe filled with a perfect fluid, in a gauge dependent 
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form/ and for all scalar, vector and tensor perturbations, neglecting anisotropic stress. 
We also present the equations in some commonly used gauges, present the evolution 
equation for a scalar field - the Klein-Gordon equation - for both a canonical and 
non-canonical field, and solve some of the evolution equations. 

First, in gauge dependent form, energy conservation at linear order gives 

5pi' + m{6p^ + (5Pi) = (po + Po)(3^i' - V'Ei' - vu!) , (3.1) 

where denotes the spatial Laplacian, = dkd^,'^ while momentum conservation 
gives 

= 0, (3.2) 

where we have introduced the covariant velocity perturbation as Vu = Vu + Bu. Note, 
Eq. (3.2) is the analogue of the Euler equation in an expanding background (see Chap- 
ter 4). 

The Einstein equations give the energy constraint equation 

SniiJi + H(j)i) - V^tpi + HEi') + nV^Bi = -AnGa^dpi , (3.3) 
and the momentum constraint 

+ \Su + Hct^u = -ArrGa^po + Po)Vu . (3.4) 



^By gauge dependent, we mean that the equations are presented in a form where the gauge functions 
a,/3 and 7* have not yet been specified 

^ Since we are working with a background whose spatial submanifold is Euchdean, the position of the 
Latin indices does not have a meaning. However, we preserve the position in order to keep with 
notational conventions such as the summation convention, and for ease of future generahsation. 



Vl, + nil-3cl)Vu + 



Po + Po 



+ 
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Finally, the component gives the equation 

+V'(0i - - + fiO + 2?/ (2^1' + V\Bi -E[) + 

= 8nGa\Po + SPi)6'j . (3.5) 

Simplifying now to scalar perturbations, the energy-momentum conservation equations 
then become 

Spi' + 3H{Spi + 5Pi) = (po + Po)(3^i' - V\E[ - vi)) , (3.6) 

6 P 

V; + n{l - 3cl)V^ + — ^ + 01 = , (3.7) 

Po + Po 

and the Einstein equations are 

SniiJi + H(t)i) - V2(^/'i + UEl) + HV^5i = -^TiGc^bp^ , (3.8) 

^\ + U4>i = -47iGa\po + Po)Vi , (3.9) 

r / 2fl"\ 

E,:'; + 2'HE,:'. + (v-i - 0i)_ ^- - - 2nB,:, + - 201 (n' - — ) 
+ v2(0i - - ^1 + 5;) + 2^1" + 27/ (^2v^i' + v'(5i - + 0/) } 

= 8nGa\Po + 6Pi)6'j. (3.10) 

We can obtain from Eq. (3.10) two scalar equations. Firstly, by applying the operator 
did^ , using the method outlined in Ref. [109], we obtain 

/2a" \ 

^Pi" + 2n^pi' + Ucpi + ( — - ^ J = 47rGa25Pi , (3.11) 
and then, taking the trace of Eq. (3.10) and using Eq. (3.11), we obtain 

{Bi - E[y + 2-H(5i - Ei) + 01 - V-i = , (3.12) 
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which can be recast in terms of the shear (ai = E[ — Bi) as 

a; + 27/(71 - 01 + = . (3.13) 

3.1.1. Uniform Density Gauge 

We now work in the uniform density gauge which, as we saw in the previous chapter, 
is specified by 6pi = 0, and consider only scalar perturbations. Eq. (3.6) evaluated in 
this gauge is 

Snm^p = (po + Po) [svi' - V\E^f! - ^p)] . (3.14) 

Noting that ipisp = — Ci is the gauge invariant curvature perturbation on uniform density 
hyper surf aces, introduced in Eq. (2.126), and 6Pisp = 5-Pnadi, this then becomes 



(Po"TPo)' 



3C( = , „, ^Pnadl - ^^E[,^ - ^p) . (3.15) 



Recall, as mentioned in section 2.4.3, that the uniform density gauge condition does not 
specify the spatial gauge function, /3i. The transformations of Ei and Vi do not depend 
upon the temporal gauge function, so Ei^p = E, and similarly for Vi. However, the 
combination Ei —vi is a gauge invariant variable, the three- velocity in the longitudinal 
gauge, and so Eq. (3.15) can be written solely in terms of gauge invariant variables as 

3Ci = -7 ^^^nadl - V^Vu ■ (3.16) 

We then recover the known result that, on large scales, the evolution of the curvature 
perturbation is proportional to the non-adiabatic pressure perturbation: 

Ci = - . '^pJ Pn.ai. (3.17) 
[po + Po) 

Thus, the curvature perturbation on uniform density hypersurfaces is conserved on large 
scales for adiabatic perturbations such as those from a single fluid or a single scalar 
field [152]. 

Let us now consider whether a general scalar field can support non-adiabatic pertur- 
bations. We can calculate the non-adiabatic pressure perturbations for a general scalar 
field with Lagrangian £ = C{X,ip), Eq. (2.41). Noting that the pressure and energy 
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density are both functions of X and (p, and so their perturbations can be expanded in 
a series as 

dP dP 

and 

Sp, = ^S^, + ^SX^, (3.19) 

the non-adiabatic pressure perturbation can then be obtained by substituting Eqs. 
(3.18) and (3.19) into Eq. (2.149), giving 

SPn.di = P,^ (^j^ - cj^ + P,x - cl^ SX, . (3.20) 

In such a system and on large scales, the relationship between Sipi and SXi is [37] 

SX, = if^Sif^ , (3.21) 

which, when substituted into Eq. (3.20), along with the definition for the adiabatic 
sound speed, c^, yields 

5Pnadl = 0. (3.22) 

Therefore, the curvature perturbation on uniform density hypersurfaces is conserved on 
large scales in an expanding universe not only for a canonical scalar field, but for any 
scalar field with Lagrangian (2.41).^ This result is in agreement with Ref. [83]. 

3.1.2. Uniform Curvature Gauge 

Now, returning to the full equations, we neglect tensor perturbations and work in 
uniform curvature gauge, where Ei = ipi = Fl = 0. The energy conservation equation 
then becomes^ 

5pi' + 3H(5pi + 5Pi) + (po + Po)vu: = , (3.23) 
while momentum conservation gives 

SPi 



v;^ + n{i-3cl)v,,+ 



Po + Po 



+ 



(3.24) 



^Note that, in a contracting universe this is not true (see, e.g., Ref. [73]). 
^We omit the tildes and gauge subscripts for brevity. 
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Let us now consider the dynamics of scalar, linear perturbations in the uniform 
curvature gauge. The energy-momentum conservation equations then reduce to 

5pi' + m{6pi + 6Pi) = (po + P^)V\i , (3.25) 
VI + H{1 - ?,cl)V, + + 01 = , (3.26) 

and the Einstein equations give 

m^(t)i + HV'^Bi = -AnGa'^dpi , (3.27) 
n(j)i = -AnGa^po + Po)Vi , (3.28) 

-Hcpi + 01 - = A-kGo^SPi , (3.29) 

B[ + 2'HBi + 01 = . (3.30) 

We will now solve this set of equations for the energy density perturbation of a perfect 
fluid. Let us first rewrite the equations in terms of the 'new' velocity perturbation 

Vi = (po + P^) (^1 + 5i) ■ (3.31) 

This enables us to write Eq. (3.28), using Eq. (1.23), the background Friedmann equa- 
tion, as 

0i = -^-Vi, (3.32) 
2po 

and using this, Eq. (3.27) becomes 

V^i?i = ^— Vi-^-5pi. (3.33) 
z Po I pQ 

Then the evolution equations, Eqs. (3.25) and (3.26) are 

5p'i + (3 + w) 6pi + m6Pi - k^Vi - (1 + w) Vi = , (3.34) 
Vi' + ^ (5 - 3w) Vi + <5Pi = , (3.35) 

where w = Pq/Po defines the background equation of state, ^ and we are working in 



'We do not demand that w be constant here. 
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Fourier space, k being the comoving wavenumber. Equations (3.34) and (3.35) make 
up a system of coupled, linear, ordinary differential equations. Given an equation of 
state and initial conditions, this system can be solved immediately (for a given k), 
numerically. One can also obtain a qualitative solution by considering a system of two 
equations such as this one. However, if one wants to solve the system quantitatively, 
and analytically, it is easier to rewrite the system as a single second order differential 
equation, which we do in the following. We solve Eq. (3.34) for Vi and get 

Vi = ^ (^5p[ + (3 + w) 5pi + mSP^ , (3.36) 

where T = fc^ + {l + w) . 

After some further algebraic manipulations of Eq. (3.36), and using Eq. (3.35), we 
arrive at the desired evolution equation 



5p'i +{7H-%\ 6p[ + (3 + w) 



^' 2:' 1 ^ ; 



6pi 

(5Pi = . 
(3.37) 



Equation (3.37) is a linear differential equation, of second order in (conformal) time. 
It is valid on all scales, for a single fluid with any (time dependent) equation of state. 
Furthermore, it assumes nothing more than a perfect fluid and hence allows for non-zero 
non-adiabatic pressure perturbations. 

Having derived a general governing equation (3.37) valid in any epoch, we now restrict 
our analysis to radiation domination, where the background equation of state parameter 
is w = I and the adiabatic sound speed is Cg = |. We recall from Eq. (2.150) that the 
first order pressure perturbation can be expanded as 



where 5Pnadi is the non-adiabatic pressure perturbation, which is proportional to the 
perturbation in the entropy, and is defined in Eq. (2.152). Then, the general governing 
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equation, Eq. (3.37), becomes, during radiation domination, where Ti oc l/r] 
3H^ 



5p'[ + 2 + 



A;2 + 6^2 



5p\ 



3'H5Pnadl + 



+ 36 



6pi 

^Pnadl =0. 

(3.38) 



For the case of zero non-adiabatic pressure perturbations the second hne in Eq. (3.38) 
vanishes, and the resulting equation can be solved directly using the Frobenius method, 
to give [lU] 



6p,ik,r]) = Ciik)r]-' 



cos I —T=] kri — 2 sin ( —'jL] 
2 cos ( V3 + sin | | kri 



(3.39) 



where Ci and C2 are functions of the wave- vector, k, and k is the wavenumber, k = \k\. 
For small kr], the trigonometric functions can be expanded in power series giving, to 
leading order, the approximation 



6pi{k, 7]) ~ A{k)k7]~^ + B{k)7]~ 



(3.40) 



for some A{k) and B{k), determined by the initial conditions. 

In order to solve for a non vanishing non-adiabatic pressure, we make the ansatz that 
the non-adiabatic pressure grows as the decaying branch of the density perturbation in 
Eq. (3.34), i.e., 

5P^Mk,v) = D{k)k^V'' ■ (3.41) 

This assumption is well motivated, since we would expect the non-adiabatic pressure 
to decay faster than the energy density, and in fact observations are consistent with 
a very small entropy perturbation today. Furthermore, if one is considering a relative 
entropy perturbation between more than one fluid or field, as time increases the system 



3.1: First Order 



59 



will equilibriate and one species will dominate.^ This gives the solution 
6piik,r]) = Ci{k)r]-'' 



cos I kri — 2 sin \ -^}=\ \/?> 



+ C2{k)r]~ 



2 cos ( ^ + (^^^) 



3D{k)k^ri~\ (3.42) 



3.1.3. Longitudinal Gauge 



In this section we consider scalar, linear perturbations in the longitudinal gauge in order 
to make a connection with the literature. In the longitudinal gauge, Ei = Bi = 0, and 
ipi = "^1, 01 = $1, and so Eq. (3.12) tells us that, in the absence of anisotropic stress, 
$1 = \I'i. Taking this into account, the energy-momentum conservation equations 
become 



6p^' + 3-H(5pi + 5Pi) = (po + Po)(3<f i + VV^ 



v[ + n{l-3cl)vi + 



6Pi 



Po + Po 



$1 = 



(3.43) 
(3.44) 



and the Einstein equations are 



3n{^[ + n^i) - V^^i = -AnGa'^Spi 



2a" 



AnGa^dPi 



(3.45) 
(3.46) 

(3.47) 



Assuming adiabatic perturbations, in which case SPi = cjSpi and w = cf, we can 
combine Eqs. (3.45) and (3.47), using Eqs. (1.23) and (1.25) to give 



$1 + 3H(i + <)$; + <r<I>i = 



(3.48) 



which, in radiation domination (where we recall from section 1.1 that = 1/3 and 
Ti oc ?7~^), then becomes 

r;$i + 4$; + \r]P<^i = . (3.49) 
o 



^Of course, given a specific model of the early universe this relative entropy perturbation can be 
calculated. 
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This equation can then be solved (see, e.g., Ref. [l u] for the general method) to give 



$l(fc,77) 



cos 



+ 



{kriY 



cos I + sin I — 7^ I kri 



(3.50) 



3.1.4. Scalar Field Evolution and Sound Speeds 

In this section we consider the dynamics of a scalar field. Many oscillating systems can 
be described by a wave equation, that is a second order evolution equation of the form 



-X- V^X + F(X,X) = 



(3.51) 



where X is, for example, the velocity potential, -F(X, X) is the damping term, and Cp^ 
is the phase speed, i.e. the speed with which perturbations travel through the system 
[79, 141]. The situation for a scalar field is not dissimilar, and the evolution equation 
in this case is the Klein-Gordon equation, where X = if is the scalar field. The Klein- 
Gordon equation is obtained by substituting the expressions for the energy density and 
the pressure into the energy conservation equation. The background equation is given 
in Section 1.2 and here we focus on the perturbations. 

There is some confusion in the literature on the meaning of adiabatic sound speed 
and phase speed. Although this seems not to affect the results of previous works, and 
the adiabatic sound speed is often simply used as a convenient shorthand for Po/po, 
as defined in Eq. (2.151), it is often confusingly used to denote the phase speed. The 
adiabatic sound speed defined above describes the response of the pressure to a change 
in density at constant entropy, and is the speed with which pressure perturbations 
travels through a classical fluid. A more intuitive introduction of the adiabatic sound 
speed is using the compressibility [79], 



/3 = 



1 dp 
pdP 



(3.52) 



which describes the change in density due to a change in pressure while keeping the 
entropy constant. 
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A collection of scalar fields can also be described as a fluid, but the analogy is not 
exact. Whereas in the fluid case the phase speed Cpjj and the adiabatic sound speed 
Cg are equal, this is not true in the scalar field case and the speed with which per- 
turbations travel is given only by c^^. The phase speed, c^^ can be read off from the 
perturbed Klein-Gordon equation governing the evolution of the scalar field. This is 
just a damped wave equation, like Eq. (3.51), and so by comparing the coefficients of the 
second order temporal, and second order spatial derivatives, we obtain the phase speed. 

We first focus on the canonical scalar field with Lagrangian (2.30); the background 
equation is given in section 1.2. The equation for linear perturbations is, in the uniform 
curvature gauge, 

Vi' + 2?/5(^; + 2a^U,^(f), - ip'oV^B^ - V^6^, - + a^U^^^6^, = . (3.53) 

We can express this in closed form, that is containing only matter perturbations, by 
replacing the metric perturbations using the appropriate field equations. Doing so gives 
(e.g. Ref. [106]) 



+ a' 



U 



8tiG /oSttG,,, , 



6(f^ = 



(3.54) 



Similarly, we can calculate the Klein-Gordon equation for a scalar field with the 
non-canonical Lagrangian, (2.41). We obtain 



c: 



ph 



3H 



c: 



ph 



5ip + 



+ C2 



5ip = 0, 



(3.55) 



where the coefficients Ci and C2, both functions of ip and X, are 
^2 



c: 



P,x 
1 



H [^oP,x^ + 'fo^P,xxJ 

P,x 



(3.56) 
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G 



P,x 

P,x^ + V^o P,xx^ (5v3o p,xx + V^o + 2p,x j 



X 



3i/(/3op,x - ^^{3H^ + 2H) + cl^{^o''p,xx^ + P,^) 



H 



H [</3o^P,Xw " P,w + S-ffv^oP.Xip] • 

P,x 



(3.57) 



The constants in this equation of motion can be interpreted physically: Gi is an addi- 
tional damping term, and G2 affects the frequency of the oscillations. The phase speed 
is then read off as 



P,x 



p^x + 2Xop,xx 

which, using Eq. (2.37) reduces to [■;(>] 

2 -Po,x 



(3.58) 



-'ph 



(3.59) 



The adiabatic sound speed for the Lagrangian (2.41) is given by 

2 P,x0o + P,v 
c = 

P,xVo - P,v + P,xx</^o Vo + P,x<^V^o^ ' 
Finally, the adiabatic sound speed for the canonical Lagrangian is 

2t/.„ 



(3.60) 



1 + 



(3.61) 



and becomes 



■1 in slow-roll. The phase speed for a canonical scalar field is 



1, as can be read off from Eq. (3.53). 



Thus we note that, although in classical fluid systems the adiabatic sound speed and 
phase speed are the same, they are different in the scalar field models studied here, 
with only the phase speed describing the speed with which perturbations travel. We 
emphasise that only the definition for the adiabatic sound speed, Eq. (2.151), enters the 
definition of the pressure perturbation. Similarly, in the adiabatic case when 5Pnadi = 0, 
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Eq. (2.152) reduces to 5Pi = 0^5 pi. A gain, this convenient relation between the pressure 
and the energy density perturbation is only correct when using the adiabatic sound 
speed, as defined in Eq. (2.151). 

3.1.5. Combined Dark Energy and Dark Matter System 

In this section, as an application of gauge choice and the linear equations, we consider 
perturbations of a combined dark matter and dark energy system. The importance of 
considering dark energy perturbations in such a system is still under discussion. Re- 
cently, in Ref. [r2()], it was shown that ignoring the dark energy perturbations results 
in gauge dependent predictions. The authors showed this by first taking the governing 
equations for the dark matter /dark energy system calculating, using a numerical sim- 
ulation, the matter power spectrum in three gauges: the comoving gauge, the uniform 
curvature gauge and the uniform expansion gauge. As expected, the same result was 
obtained for each gauge. They then ignored perturbations to the dark energy compo- 
nent, and repeated the analysis, this time obtaining different results for each of the 
three gauges. They then concluded that it was crucial to include the dark energy per- 
turbations when analysing such a system. Thus one cannot use the familiar evolution 
equation for the dark matter density contrast, 6c = 6pc/pc, 



which is obtained by ignoring dark energy perturbations once the equations have been 
put into a gauge comoving with the dark matter. 

In this section we consider a similar question, and ask whether we can come to the 
same result without using numerical simulations, but instead by using the formalism of 
cosmological perturbation theory. The section contains work published in Ref. [V>]. 

We consider here the linear governing equations for scalars only in coordinate time 
and without fixing a gauge (dropping the subscript '1' in this section for ease of presen- 
tation). Though some of the equations from previous sections are replicated here, since 
we are using a different time coordinate in previous sections, we give all equations used 
for completeness. We assume that the dark matter and dark energy are non-interacting, 
and so energy momentum conservation for each fluid is given by 



4 + 2H5c - 47rGpc^c = , 



(3.62) 



0. 



(3.63) 
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Note that this assumption is for brevity, and to allow us to deal with more manageable 
equations. The results highlighted in this work will still hold in the case of interacting 
fluids for which the overall energy-momentum tensor is covariantly conserved, but the 
components obey 

V.Tl"^^, = Q^o.)u , (3.64) 

where Q[a)iy is the energy-momentum transfer to the ath fluid [70, 112]. 

Then, we obtain an evolution equation for each fluid from the energy (temporal) 
component of Eq. (3.63), 

5p^ + 3H{5pa + 5Po) + (po + P^)^{Va + cr) = 3(p, + P„)^ , (3.65) 

where the covariant velocity potential of each fluid is defined as = a{va + B). 

Considering the dark matter fluid and the dark energy scalar field, respectively, 
Eq. (3.65) then gives 

4 = 3^-— (a + K), (3.66) 



a2 



5^ + ?,H5^ + {u,^^ -^y^ = ^ {?,^ + 2f/,^0 , (3.67) 

where we have used the expressions for the energy density perturbation and pressure 
perturbation of a scalar field given in Section 2.2, and note that the dark matter is 
pressureless. We have also used the background evolution equation for the dark energy 
scalar field 

00 + 3if¥3o + f/,^ = , (3.68) 
the background conservation equation for the dark matter, 

Pc + 3Hpc = , (3.69) 

and the expression 

V; = -— . (3.70) 
There is also a momentum conservation equation, coming from the spatial component 
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of Eq. (3.63), corresponding to each fluid 

K - ^HclV^ + + ^ = , (3.71) 

Pa + Pa 

where 4 = Pa/ Pa- 

The Einstein field equations give (from the previous section or, e.g., Ref. [Ill]) 

?,HN) + H(P) -N) + Ha)+AnG6p = Q, (3.72) 

ij + H<p + A'KG{pQ + Po)V = (3.73) 
(T + i7(T-0 + V^ = O, (3.74) 
ti^ + ZHi; + H<j)+ {3H^ + 2H)(j) - 47rG6P = , (3.75) 

where the total matter quantities are defined as the sum of the quantity for each fluid 
or field, i.e. 

Sp = Sp^ + Spc, (3.76) 
6P = 6P^, (3.77) 
(p + P)V = {p^ + P^)V^ + p,V,, (3.78) 

and we have used the fact that the dark matter is pressureless, i.e. Pc = = SPc- 
Introducing a new variable Z, both for notational convenience and to assist with the 
following calculations, defined as 

Z = 3ii; + H^)-—a, (3.79) 
we can rewrite Eq. (3.75), using Eqs. (3.72) and (3.74), as 

Z + 2HZ + [3H + — j = 47iG{Sp + 3SP) . (3.80) 



Then, from Eq. (3.66), 

Z = 5, + 3H(j)+\v,, (3.81) 
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and from Eq. (3.71) for the dark matter fluid (for which Cg 



0), 



= -K. 



(3.82) 



Differentiating Eq. (3.81) gives 



Z = 6, + 3{H(i))- + 




) 



(3.83) 



Substituting this into Eq. (3.80) gives 



+ V;(6(if + H"^) + 167rGv3o^) + 3i/K . 



(3.84) 



The evolution equation for the field is then obtained solely in terms of matter pertur- 
bations from Eq. (3.67) by using Eqs. (3.82) and (3.81): 



It is worth restating that we have derived these equations in a general form without 
fixing a gauge. If we set the dark matter velocity to zero, i.e. = 0, then they reduce to 
those presented in, for example, Refs. [34, 70] (for the case of a zero energy-momentum 
transfer) . 

We now want to consider fixing the gauge freedoms in this set of governing equations. 
In order to do so, we need to consider the gauge transformations of the variables under 
the transformations 




(3.85) 



X* = X* - /3/ , 



t = t — a , 



(3.86) 
(3.87) 



where the generating vector is then defined as 



e = («,AO 



(3.88) 



3.1: First Order 



67 



Note that, since the choice of time coordinate is different to that used in the previous 
chapter, the following gauge transformations will differ from those presented in, e.g.. 
Section 2.3.4 due to the different definition of Then, scalar quantities such as the 
field perturbation transform as 

6ip = 6ip + ipa , (3.89) 

the density contrast transforms as 

4 = <5c + , (3.90) 
and the components of the velocity potential as 

Va = Va — CLf3 . (3.91) 

Furthermore, by considering the transformation behaviour of the metric tensor we ob- 
tain the following transformation rules for the scalar metric perturbations 

= + d, (3.92) 
^ = ^-Ha, (3.93) 

B = B --a + a$, (3.94) 
a 

E = E + /3, (3.95) 
and so the scalar shear transforms as 

a = a + a. (3.96) 

Finally Eq. (3.91) gives the transformation behaviour of the components of the covariant 
velocity potential 

Va = Va-a. (3.97) 

Turning now to the case at hand, choosing the gauge in which the perturbation in 
the dark energy field is zero, 6ip = 0, fixes a as 
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Since none of the gauge transformations of the quantities involved in the governing 
equations depend upon the threading /3, we do not need to consider fixing it exphcitly 
here. (Of course, one can rigorously fix (3 by choosing a suitable gauge condition, as 
outlined in Section 2.4.) The governing equations in this gauge are then 

4 + 2H6, - AnGpX = V,{6{H + H^) + IQ-kG^"^) + ^HV^ , (3.99) 

K - — - = -Sc , (3.100) 

where the hat denotes that the variables are evaluated in the uniform field fluctuation 
gauge. That is, in this gauge, 6c and Vc are gauge invariant variables defined as 

6c = 6c-^6ip, i?=K + ^. (3.101) 
PcV^o V^o 

Alternatively, choosing a gauge comoving with the dark matter, in which Vc = fixes 
the generating vector as 

a = Vc, (3.102) 

and reduces the governing equations to 

(fc + 2Hic - 47rGpc4 = 87rG'(2v3o(5y? - U^^Sip) , (3.103) 
5^ + 3HS^ + [U,^^ -—y^ = v3o4 , (3.104) 

where the bar denotes variables in the comoving gauge and we have 

4 = 5c + -K, 6'ip = 6ip + ipoVc. (3.105) 

Pc 

By studying the above systems of equations, it is evident that choosing the dark 
energy field perturbation to be zero is a well defined choice of gauge, reducing the 
governing equations to Eqs. (3.99) and (3.100). Then, having done so, we are no 
longer allowed the freedom to make another choice of gauge. Alternatively, choosing 
a gauge comoving with the dark matter uses up the gauge freedom, and so we are 
not permitted to neglect the perturbation in the dark energy field. In fact, doing so 
will result in erroneous gauge dependent results. It is clearest to see why this is the 
case by considering the set of governing equations. By making our choice of gauge we 
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are left with a set of equations which is gauge invariant: that is, performing a gauge 
transformation will leave the set of equations unchanged. However, by neglecting the 
perturbation in the dark energy field after having chosen the gauge comoving with the 
dark matter amounts to setting the right hand side of Eq. (3.103) to zero. This resulting 
equation will, in general, then no longer be gauge invariant. 

Thus, we conclude that the dark energy perturbation must be considered in a system 
containing a mixture of dark matter and dark energy. Our result is consistent with that 
of Ref. [126], though we have shown this by simply using the formalism of cosmological 
perturbation theory instead of relying on more involved numerical calculations. 



We can extend the governing equations presented in the previous section to beyond 
linear order by simply not truncating the expansion of each variable after the first 
term. Doing so, we obtain equations with similar structure to those at linear order, 
however with new couplings between different type of perturbation. In fact, these 
couplings will turn out to be the reason for the qualitative difference between the 
linear and higher order theories. In this section, we will present the full second order 
equations for scalar, vector and tensor perturbations in a gauge dependent format. 
These equations have been derived previously, for example in Ref. [123], though are 
derived in a slightly different way, using the ADM split, and so are written in terms of 
the extrinsic curvature. The energy conservation equations have also been derived in 



Conservation of energy momentum again gives us two equations: an energy conser- 
vation equation 
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Ref. [112]. 



5p^ + 3H(5p2 + 5P2) + (po + Po)(C^: + V2^:) + 2(5pi + 5Pi)^^\ 
+ 2(5pi + 5Pi)(Cf, + vu!) + 2(po + Po) [{Vi + v\)Vi, + 
- 2C[,f'^ + v\{C{,, + 20i_J + mv\{Vu + i^i.)! = , 



(3.106) 
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and a momentum conservation equation 

(Po + Po)V2^] ' + (po + Po) (02, + mV^,) + + 2 [{6p^ + 6P,)Vu 
+ 2{6p, + 5Pi)(0i,, + AHVi,) - 2(po + Po)' [(^i^ + ^i.)0i - ^C^X 
+ 2(po + Po) \Vu{Ci'^ + vl,) - Bu{<pi + m<Pi) + {2CuA 



0, 



(3.107) 



where we have not expanded the terms in Cij in order to keep the equations compact, 
but recall from Eq. (2.3), that 



Cij = —ijjSij + E^ij + + -hij . 
The Einstein equations give the (0 — 0) component 



(3.108) 



v'ci, - C2.,;^ + 2H(-q , + Bl, + m<P2) + 2Cij,{-cik: - 2Cf ,,) 

+ 2B{ \c(,, - C[,^; + \ [V'Bu - By/) + 2H [C^^ - 2Cu,^ - 



2C\jk^i' — C^k^ij — V'^Ciij + 2'H{C^^j — Biij) + 2Cijk/Cf' 



+ 2Cl .(P,/ - -CI , + m^) + .C,,^ + 2C,,(-Cr - Bi 
+ ^5i,,(5i/' + - Qn\m - BuB\) - B{,B,^,^ - mB\,ct>^ 



-SnCa^ 



2(po + PoWl'vik + 5p2 



(3.109) 
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and the (0 — i) component 



2 (Si,,, - Si,,,) Qcffe,^' - C^^fc^ - 2Ci,fe,, [Bl^ - S^,'^) + 2S^,,0i, 



+ 



Si,y-V2SH + 2(C;,/-Ci^' 



leyrC 



:F2^ - 0i(V^i. + Si,) + 2Cukvl + (5pi + 5Pi)Vi 



li 



(3.110) 



We present the full {i — j) component in Appendix A. 



The equations for scalar perturbations only in a gauge dependent form are then 
obtained by substituting C^,- = —ipSij + E^ij and S^ = S,^, at both first and second 
order, into the above. The energy conservation equation then becomes 

5p2' + m{5p2 + 5P2) + (po + Po) {y\E^ + V2) - 3^2') + 2(5pi + 5Pi),^i' 
+ 2(<5pi + 5Pi) (v'(Si' + vi) - 3^1') + 2(po + So) [{VI, + 4-Ht;i,,)(\/i,, + t;i,,) 
+ 3V'iV^i' + V2t;i0i - (V^iV^^S)' + S(_,^.Si,^^' + t;i,*(20i,, - 37/^1,, + V^Si,, 



= 0, 
(3.111) 



while the momentum conservation equation is 

(po + P^)Vi,^ ' + (po + So) (02 + mV2) ^ + 5P2, + 2 [Vi,,(5pi + 5Si 
+ 2(5pi + 5Si) f 01 + AHV^ - 2(po + So)' f(V^i + 5i),,0i - 2(Si,„t;i,^' - ^it;i,,) 



+ 2(po + So) V^i,,( V'(Si' + v^) - 3^1' ) - Si,,(0i' + m<Pi) + t;i,^(t;i,„ + 8?^Si 



+ 2(VSi,i,-^it;i,,) - 01 ( (Vi + Si)' + 201 + 4?/t;i j ^ - 8-H^it;i,, =0. (3.112) 
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Turning now to the Einstein equations, the energy constraint is 

3niii^2' + ncj)2) + (niB2 - e^) - + v'^i (v'(^i' - ^Si) - 

+ Si,, (n{3nB^; - 2w''e,; - 2(V'i + 0i);) - 2V'i;') + 2e^;^{xIj, - 2nB^),,^ 

+ 47/(^1 - 0i) (3^1' - V\E^ - El)) + E^;^' (mE^ + ^Ei' - El) , , 

+ V'i'(2V'(£;i' - 2^^) - 3V'i')) + V'i;(2V'£;i - 3V'i),i + 2VVi(V'£;i - 4V'i) 

- 12?^^0i' + ^ (Ei,,,£i,'^' + V'Ei,,V'Ei/' - E.^i^kEi;^" - V'^i'V'Ei') 

= -4776-02 (2(po + Po)T4>i,fe + <5P2) , (3.113) 

and the momentum constraint 

i^2,i + nh, - + 01 + v^Ei)/ + B,^,pnB, + 0i)/ 

- V'£;i - 4V'i)] ' - 01,, (87^01 + 2V'i' + V\E^ - bS) 

- B.ji;,/ + 2ij,f'E,^,j + E.j.'E,/'' - ^/.(V^^i + 40i) - V'^i^i,, 
= -iirGa^ [(po + Po) (v^2,i - 20i(Vi + B,),i - 4(^ii;,, - E.^.kVi,")) 

+ 2(5pi + 5Pi)yi,,] , (3.114) 
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while, from the trace of the i — j component, we obtain 

3?^(2^2 + 02)' + V\E2" + 2E2' + 2V^2 - B2' - 02 + 2nB2) - 302 (n^ - 2—) + 3V^2" 
+ (V'l - 0i) (l2(V'i" + 27^^'/) + 4V'(0i + (Si - + 8?{V'(Si - 

+ Ei;^' (8H(i?i' - fix),, + 2^1,,^. - 401,^. - 451^) + E^:^ (^E/,^. - 5i,,.) 

+ Tsi'^E^ (40/ - v2(£;i' - 2Si)) + v^£;i; (v2£;i,, + 20i,, - mB^^^ - 25/,) 

+ (2V'£;i,, - 4?{5i,, - 2(^1 + 0i),, - 25i;,) - 20i'(V25i + 12?/0i) 

- 201^,01; + 2vVi(v'£;i - + ^i'(3v^i' - 60i' - 8w'£;i - 2v'(£;i' + ^o) 
+ 25i; (?^(35i' - 20i) - 3^1') ^ + ^ - - v'^iV'Si) 

+ 4(£;i;^£;i;^ - ^^'V'e^) + 3(n'- 2^) (40^^ - b,,,b,;) 

= 47rGa' (35P2 + 2(po + i'o)t'i,'V^i,i) • (3.115) 

We can obtain, from the i — j component a fourth Einstein equation, as at hnear 
order, by applying the operator did^ . However, since we do not require this equation 
for the work in this thesis, as one can always use the energy-momentum conservation 
equations in place of the i — j equations, we omit that equation here. 

3.2.1. Uniform Curvature Gauge 

Wc now present the second order equations in the uniform curvature gauge, where 
-0 = £^ = Fj = 0, at both first and second orders. We also consider only scalar and 
vector perturbations in this section, i.e. we choose to neglect tensor perturbations, and 
so Cij — 0, at both first and second order. We obtain the energy conservation equation 

Sp2' + m{5p2 + 5P2) + (po + Po)v2i: + 2(5pi + 5P^),^{ + 2(5pi + 5P^)vu: 

+ 2(po + Po) \{Vi + v\)Vi^ + vii<t>^ + 2^;i0i,, + 4^1(^1, + v^)] = , (3.116) 
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and momentum conservation equation 



(Po + ^^0)^2.] + (Po + + mV2i) + + 2 [(Vi + 5P^)Vu 

+ 2(5pi + 5Pi)(0i,i + 4?^Ki) - 2(po + i"o)'(5H + 
+ 2(po + Po) Wuvlj - Bu{<t>i + 87^01) + v{{Vu,j - 



0i(V^i', + i?i, + 20i, + 4?^^;H) 



0. 



(3.117) 



The Einstein equations give us an energy constraint 



2(po + Po)V^vik + 5P2 



(3.118) 



a momentum constraint 
1 



Bu,, + + 2fi5',,0i,, + 201 (Si,/ - V'Sh) 

1. 



IGTrG 



(3.119) 



and a third equation from the {i — j) component 



\ a 



n^) + 2H{Blk + 0^) + + VV2} + 20i;0ij 



+ {Bij^k + B\jk — 2-BiA:,*j) + (-BJ J + -Bij/)(0i' + -Bi,fc) — Bik^'B^j — Bij^ B\^k 
+ B\{B,,/ - V^Si, + 47^01,,) + 20i[Pi,;' + B,^\, + 201/, + 2n{B,,: + Pi\,)] 



2a" 



(40? - SifeSi*^) - 20i[Sf - V'0i + 27^(20; + Blk)\ 



1 



-Bi [V^-Bife - -Bii,A:' + 2H{B[^ - 0i_jk)] - -{2B'1 ^kBiil' - -Bi/,fe-Si/ - 3-B{,fe-Biz,'^) 



0iPi',fe - 0i,fe0i,'} = SttG [2(po + Po)v\V,j + 5'j5P2\ . 



(3.120) 
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Considering scalar perturbations only, energy conservation gives 

bp^ + m{bp2 + + (Po + Po)V2t;2 + 2(<5pi + 5Pi),iT;i' 
2(5pi + 5Pi)V'^; + 2(po + Po) [(n'i + 4?^^i,0(^i,i + 



+ V'^;l0l + 2T;i,Vi,iJ = , 
and momentum conservation gives 

(PO + Po)^l,i] ' + (po + Po) (02 + 47^1^2) . + ?>P2,i + 2 + (^Pl) 

+ 2((5pi + (5Pi) (01 + 4?^yi) . - 2(po + Po)'(T4 + 5i) 
+ 2(po + Po) k,iV'^;i - Pi,i(0i' + m4>x) + 



(3.121) 



(/)i((T/i + Piy + 20i + 4?{^;i) 

\ / , 1 



0. 



(3.122) 



Finally, the Einstein equations for scalars only in the uniform curvature gauge are 

37^202 + V'T^Ps J(V'Pi)' + UB^pUB^^ - 20i;) + 47{0iPi 

- 12?^Vi' + ^^i.^i^i/^' = -47rG'a2 (2(po + Po)V^i, + ^^2) , (3.123) 



and 



n4>2,i + B.^.pUB^ + 00/ - 01,, (^87^01 - V^^i^ 
= -47rGa' [(po + Po) (^2,i - 201(14 + Pi),,) + 2(5pi + 5P^)V^, 



3'H(p2 + V\2'HB2 - B2' - 02) - 01 (4V='(0i + Pi') + mV^B^ 
- 20i'(V'Pi + 12?{0i) - 20i,,0i; + 27{Pi;(3Pi' - 20i) 



+ oiPhijB. 



(3.124) 



= 477^0" (35P2 + 2(po + Po)vi%,^ . (3.125) 
The Klein-Gordon equation can be obtained at second order by using the same tech- 
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nique as at first order: comparing the energy-momentum tensor for tlie scalar field to 
that of a perfect fluid and using energy conservation. The equation for a canonical 
scalar field is 



2Vi(V^5i + 



4Vi .i^i,' = , 



(3.126) 



where we are yet to use the field equations to remove the metric perturbations. See 
Refs.[65, 67, 68, KK ] for the closed form of the Klein-Gordon equation at second order 
and for detailed work on the second order Klein-Gordon equation. 



3.2.2. Poisson Gauge 

In this section we present the second order equations in the Poisson gauge, in order to 
connect with the literature which often uses this gauge (for example, Ref. [3] presents 
the Einstein equations with scalar field matter in this gauge). The gauge is defined by 
E = = B, and then = $ and ip = In the absence of anisotropic stress, as is the 
case for this work, \I'i = $1 (though note that this does not hold true for the second 
order variables $2 and \E'2)- Note also that, in this gauge, V = v. 
Energy conservation then becomes 

Sp2 + m{5p2 + 5P2) + (po + Po) {^^V2 - 3^'2) + 2((5pi + 5Pi),^vi' 
+ 2(5pi + 5Pi) f V V - 3$;) + 2(po + Po) b{v[^i + mv^^i)v^: 



0. 



(3.127) 



while the momentum conservation equation is 

(PO + P0)V2^ ' + (po + Po) f $2 + mV2 ) + 5P2,^ + 2 Vi^i{5pi + 5Pi] 



+ 2(5pi + 5Pi) ($1 + mvr^ ^ - 6(po + Po)'$it;i,. 

+ 2(po + Po) \vi,, - 3$;) + t;i>i,,, - $1 (v[ + 2<l>i + AUvi 



(3.128) 
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Then, the Einstein equations (where we here do not decompose the velocity into a 
scalar and divergenceless vector) are 

3H{% + H^2) - v2^2 - 3$;$; - 3<i>i;<i>i,i - sv^^i^i - uu^^i 

= -A-KGa" (2(po + PoWv^k + 5P2) , (3.129) 



^2 . + ^^2,i + 4($i,,<i>)' - <^iA^u^i + 2$;) - 4$; ,$1 



(PO + PQ){V2^ - Q^lVvd + 2(5pi + 5Pi)Vii 



(3.130) 



and 



1 / 2a" 

% + H(2^2 + $2)' + 3 V'(<l>2 - ^2) + ( — - ) $2 



V a 



2a" 



2<i>i;<i>i,, - 87/$!$; - -v^$i$i - 3($;)2 



47rGa' 5P2 + 77(po + PoWvu 



(3.131) 



For completeness, we present the fourth field equation, obtained by applying the oper- 
ator did^ to the i — i component of the Einstein equations, Eq. (A.l): 



2a" 



hA^2 = ^7iGa^6P2 + 8nGa\po + Po)V-^did^ {vi'vi^ 



V-^<l 2$i + 4$i \<l>i - V 



. 3 



2a" 



(3.132) 



where V ^ is the inverse Laplacian operator. Finally, combining Eqs. (3.131) and 
(3.132), we obtain 



V'(*2 - $2) = 247rGa2(po + Pq) [v\vu - V'^ {8^8^ {v'v,))^ + 12^l(n^ 

- i8H$i$; - 3V"2|2<i>i/$i,'=' + 4<i>i;j$i/ + <i>;(4vVi - 3v2$; + 2hv2$i)| 



2a;^ 

a 



- 6$i,i<i>i; + $1$" + ($'1)2 + 2<i>2 (n^ 



2a" 



(3.133) 



which is the second order analogue of the equation which, at first order, tells us that 
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the two Newtonian potentials are identical in the absence of anisotropic stress. 

3.3. Discussion 

In this chapter we have presented the dynamical equations for general scalar, vector 
and tensor perturbations of a flat FRW spacetime at both linear and second order. The 
case of linear perturbations is relatively simple since the different types of perturbation 
decouple from one another. However, as one moves beyond linear order, this is no 
longer the case and so things necessarily become more involved. As we have seen, 
at second order the energy-momentum conservation equations, for example, not only 
depend upon the true second order perturbations, but also involve terms quadratic 
in first order perturbations. So, while in this chapter we have managed to solve the 
linear equations analytically in, for example, the uniform curvature gauge, doing this at 
second order would be far more complicated. But with this complication comes great 
reward, since this coupling between lower order perturbations provides a source which 
can result in qualitatively new results, and thus new observational phenomena. In the 
next chapter we will discuss one such example: vorticity generation at second order in 
perturbation theory. 



4. Vorticity 



Vorticity is a common phenomenon in situations involving fluids in the 'real world' (see 
e.g. Refs. [2, 79]). There has also been some interest recently in studying vorticity 
in astrophysical scenarios, including the inter galactic medium [14, 155], but relatively 
little attention has been paid to the role that vorticity plays in cosmology and the early 
universe. 

In this chapter we will consider vorticity in early universe cosmology. Starting with a 
summary of vorticity in classical fluid dynamics as motivation, we show that vorticity 
is generated by gradients in energy density and entropy. We then consider vorticity in 
cosmology, for which we need to use general relativity and cosmological perturbation 
theory. At linear order, there is no source term present in the evolution equation, and 
any vorticity present in the early universe will decay with the universe expansion. At 
second order, however, vorticity is induced by linear order perturbations. As mentioned 
in Chapter 1, while at first order perturbations of different types decouple, this is no 
longer true at higher orders. Recent work in the area of second order gravitational waves 
has exploited this fact [11, 13-15, 20, 114, 117, 125, 132, 148], as have recent studies 
of induced vector perturbations [97, 98]. Though Ref. [98] assumed the restrictive 
condition of adiabaticity which therefore could not source vorticity at any order, we 
show that, in analogy with the classical case, vorticity is sourced at second order in 
perturbation theory by a term quadratic in energy density and entropy (or non-adiabatic 
pressure) perturbations. Finally, we present a first estimate of the magnitude and scale 
dependence of this induced vorticity, using the expression for the energy density derived 
in Section 3.1.2 as an input power spectrum along with a sensible ansatz for that of 
the non-adiabatic pressure perturbation. We close the chapter with a discussion of the 
results, and highlight some potential observational consequences. The results in the 
chapter have been published in Refs. [38-40]. 
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4.1. Introduction 

In classical fluid dynamics the evolution of an inviscid fluid in the absence of body 
forces is governed by the Euler, or momentum, equation [79] 

Q -I 

-^ + {v-V)v = --VP, (4.1) 

where v is the velocity vector, p the energy density and P the pressure of the fluid. 
The vorticity, cj, is a vector field and is defined as 

u} = Vxv, (4.2) 

and can be thought of as the circulation per unit area at a point in the fluid flow.^ 
An evolution equation for the vorticity can be obtained by taking the curl of Eq. (4.1), 
which gives 



V X (v X w) + — Vp X VP . (4.3) 



- = Vx(.xc.) + - 

The second term on the right hand side of Eq. (4.3), often called the baroclinic term in 
the literature, then acts as a source for the vorticity. Evidently, this term vanishes if 
lines of constant energy and pressure are parallel, or if the energy density or pressure 
are constant. A special class of fluid for which the former is true is a barotropic fluid, 
defined such that the equation of state is a function of the energy density only, i.e. 
P = P{p), and so (l/p2)Vp x VP = 0. 

For a barotropic fluid, the vorticity evolution equation, Eq. (4.3), can then be written, 
by using vector calculus identities, as 

— = — + {vV)u = {u-V)v-u{V-v), (4.4) 



^The circulation, F, is defined as 

T = <j> v dl, 

where C is the boundary of the surface 5". Then, using Stokes theorem, this becomes 

r = j j {V xv)-dS, 

which gives the result that the vorticity is the circulation per unit area at a point in the fluid flow. 
□ 
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where D/Dt denotes the convective, or material derivative, which is commonly used in 
fluid dynamics. From Eq. (4.4) it is clear that the vorticity vector has no source, in 
this case, and so c*; = is a solution. 

Allowing for a more general perfect fluid with an equation of state depending not 
only on the energy density, but of the form P = P{S, p) will mean that, in general, the 
baroclinic term is no longer vanishing, and so acts as a source for the evolution of the 
vorticity. This is Crocco's theorem [13] which states that vorticity generation is sourced 
by gradients in entropy in classical fluid dynamics. 

4.2. Vorticity in Cosmology 

In General Relativity, the vorticity tensor is deflned as the projected anti symmetrised 
covariant derivative of the fluid four velocity, that is [75] 

where the projection tensor V/j^u into the instantaneous fluid rest space is given by 

'PiM,. = QiJiu + Uf,Uy. (4.6) 

Note that, in analogy with the classical case, it is possible to deflne a vorticity vector 
as uj^ = ^e^y-yOJ^'^ , where e^^,^ = u^Ss^y^ is the covariant permutation tensor in the fluid 
rest space (see Refs. [62, 98]). However, since this is less general we choose to work with 
the vorticity tensor when deriving the equations, and only switch to using the vorticity 
vector when solving the equation in Section 4.4. 

The vorticity tensor can then be decomposed in the usual way, up to second order in 
perturbation theory, as Uij = couj + \uj2ij- Working in the uniform curvature gauge, and 
considering only scalar and vector perturbations, we can obtain the components of the 
vorticity tensor by substituting the expressions for the fluid four velocity, Eq. (2.23), 
along with the metric tensor into Eq. (4.5). At flrst order this gives us 



(4.7) 
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and at second order 



The first order vorticity is gauge invariant. In order to see this we recall, from Eq. (2.94), 
that Vii transforms under a gauge transformation as 



Vii = Vii - ai, 



(4.9) 



so the first order vorticity transforms as 



(4.10) 



and is therefore gauge invariant. 



4.3. Vorticity Evolution 

In order to obtain an evolution equation for the vorticity at first order, we take the 
time derivative of Eq. (4.7) to get 

=«%M] + «K[M]- (4-11) 

Noting that, from Eq. (3.24), 

V;'f,,.]+?/(l-3aVi[„.] + 

which gives 

Vli^,] = -ni - 3c^)Vi[,,] = -^mi - 3cj)uu, , (4.13) 
and so, from Eq. (4.11), 

u[,^-3ncluu, = 0. (4.14) 

This reproduces the well known result that, during radiation domination, {uujU^i \ oc a^^ 
in the absence of an anisotropic stress term [7-^)] . 

At second order things get somewhat more complicated. We now take the time 



+ 



PO + -Po . Mi] 



0, 



(4.12) 
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derivative of Eq. (4.8), to give 



^2ij = (^'^2[i,j] + Oy2[i,j] + 2a' V[y^ij] + {Vi + Bi).^ - 

+ 2a [v/f,Vi,] + ^/[.V^/^.] + (Vi + 5i)^.] + (^1 + ^O;-] 



(4.15) 



Therefore we now must use the first order conservation and field equations to eliminate 
the first order metric perturbations as well as the second order conservation equations in 
order to eliminate the second order metric perturbation variables. This process involves 
simple algebra, but is very tedious and so we omit the intermediate steps and instead 
quote the result. We arrive at the evolution equation for the second order vorticity 



UJ 



2ij 



5Pi + 5pi^ ^ ^ 



Po + Po 

+ 2 {V,' - X'^) ujuj,k - 2 {Xl^ - V,'J uuk + 2 [Xl^ - V^,) u^^u 

ST-L (V^li^-Pnadlj — VijSPna.dl,i) 
1 



(4.16) 



Po + Po 



+ 



Po + Po 

where Xn is given entirely in terms of matter perturbations as 



Xu = V- 



{5pi^,-nipo + Po)Vu) 



(4.17) 



Eq. (4.16) then shows that the second order vorticity is sourced by terms quadratic in 
linear order perturbations. 

In fact, even assuming zero first order vorticity, i.e. Uuj = 0, the second order 
vorticity evolves as 



^211 ~ 'i'HcluJ2ij 



2a 



Po + Po I 



i^-Pnadlj] + 



^Pl,[j^Pnadl,i] 

Po + Po 



(4.18) 



and so we see that there is a non zero source term for the vorticity at second order 
in perturbation theory which is, in analogy with classical fluid dynamics, made up of 
gradients in entropy and density perturbations. Note that, in the absence of a non- 
adiabatic pressure perturbation, we recover the result of Ref. [9S] that there is no 
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vorticity generation. 



4.4. Solving the Vorticity Evolution Equation 

Having derived an evolution equation for the second order vorticity in the previous 
section, we now seek an analytic solution to this equation (or, more precisely, the 
power spectrum of the vorticity governed by this evolution equation), using suitable, 
realistic approximations for the input power spectra. 



4.4.1. The Vorticity Power Spectrum 

In order to keep our results conservative and our calculation as simple as possible and 
hence analytically tractable, we assume that the source term in Eq. (4.18) is dominated 
by the second term.^ Then, choosing the radiation era as our background in which 
Cg = 1/3, the evolution equation simplifies to 



9a 



- Hujij = ^Sp^y6Pn^d,i] , (4.19) 



V5 



where we note that, for the remainder of this Chapter, we omit the subscripts denoting 
the order of the perturbation in order to avoid notational ambiguities, and to keep the 
expressions as compact and clear as possible: the vorticity is a second order quantity and 
the energy density and non-adiabatic pressure perturbations are first order quantities. 
We define the right hand side of Eq. (4.19) to be the source term, 

9a(n) 

%(^,^) = ^;^^P,b-^^nad,], (4.20) 



which on defining the function f{ri) as 



9a 
16pg 



fiv) = T7rr2 ' (4-21) 



is written, for convenience, as 



S,j{x, f]) = 2/(r7)5p,y5P„ad,i] . (4.22) 

^This is a reasonable assumptions since we are working on small scales and the first term has a 
prefactor 
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Since we want to solve the evolution equation we move to Fourier space, in which the 
source term becomes the convolution integral 



Si,{k,r)) = -T^^f^ / d'Hkkj - ~kjki)6Pr,,d{k,r])5p,{k - k) , (4.23) 



m 

(27r)3/2 

where k is the wavevector, as usual. Instead of considering the vorticity tensor Uij, it is 
easier, and more natural in this case, to work, in analogy with the classical case, with 
the vorticity vector. This is defined as^ 

(Ji{x, r]) = eijkUJ^''{x, 7]) , (4.24) 

where eijk is the totally antisymmetric tensor, and we can define a source vector in an 
analogous way: 

Si(x,ri)^eijkS'''(x,ri). (4.25) 

Note that the vorticity is an axial vector (that is, it arises from the generalisation 
of the cross product), and so both Ui and Si arc pscudovcctors. Thus, under the 
transformation of the argument x — )■ —x the vector vorticity transforms as tUj —i^i, 
and similarly for Si. The source vector is then Fourier transformed as 

with the source vector in Fourier space then being split up as 

Si{k,7^)=SA{k,7i)et, (4.27) 
where Sa are the amplitudes, with A e {1, 2, 3}, and the basis vectors are 

{el el e?} = |e,(fc), e,(fe), ^ (4.28) 

In order to keep a right handed orthonormal basis under the sign reversal of fc (fc — > — fc) , 



^Though strictly this is a covector not a vector, since we are working in a flat background the two are 
equivalent up to raising or lowering of indices by the Kronccker delta. Therefore, we are slightly 
loose with terminology here, and do not differentiate between the two. 
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the basis vectors must obey 

ei(-fc) =ei(fc), (4.29) 

e,(-fc) = -e,(fc) ; (4.30) 

the basis vectors are also cyclic: 

eijke{el = C'ii. (4.31) 

Given these definitions, the evolution equation, Eq. (4.19), can then be written as 

uj'A{k,r])-nuJA{k,ri)=SA{Kr]), (4.32) 

for each basis state, A (though we omit the subscript in the next few lines, since the 
evolution equations is the same for each polarisation). The left hand side of Eq. (4.32) 
can be expressed as an exact derivative, giving 

aHMl\^S(k,r,), (4.33) 

which, in radiation domination when a = t], becomes 

Uj{k,T]) 



V 

This can then be integrated to give 



7]~^S{k,r]). (4.34) 



rv 

uj{k,r]) = 7] I fj^^S{k,fi)dfi , (4.35) 
•J no 

for some initial time rjQ. Having solved the temporal evolution of the vorticity, we now 
move on to considering the power spectrum. In analogy with the standard case for 
scalar perturbations, we define the power spectrum of the vorticity as 

{uj*{ki, 7])u{k2, T])) = ^5(fci - k2)VUk, v) , (4.36) 
where here the star denotes the complex conjugate, and k = \k\ is the wavenumber, as 
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usual. On substituting Eq. (4.35) into Eq. (4.36), we can write the correlator for the 
vorticity as 



rv rv 

{uj*{ki,r])uj{k2,7])) = r]^ r^fW rj^^ {S*{ki,fii)S{k2,f]2)) driidfi2 



(4.37) 



and so, in order to obtain the vorticity power spectrum, we must calculate the correlator 
of the source term. Thus, we need to consider how the Fourier amplitudes, SA{k,T]), 
behave under complex conjugation. From Eqs. (4.26) and (4.27) we can write 



1 



(27r)3/2 

whose complex conjugate is then 
1 



Si{k, 7]) ei{k) + S2{k, r])ei{k) + 53(fc, ifjh 



ikx 



(4.38) 



(27r)3/2 



d'k 



Sl{k,r]) ei{k) + Sl{k,vi)ei{k) +S^{k,r])ki 



-ikx 



(4.39) 



Alternatively, under the change k —k, on which — )■ —Si, as mentioned above, 
Eq. (4.38) becomes 



Si{x,r]) 



1 



(27r)3/2 



d^k 



- Si{-k, 7]) ei{-k) - 52(-fc, r])ei{-k) 
- Ssi-k^T]) ■ {-h 



-ikx 



(4.40) 



Comparing Eqs. (4.39) and (4.40), we can then read off the relationship between the 
Fourier amplitudes and their conjugates: 



Slik,r]) = -S^i-k,r^) 
S;{k,r]) = S2{-k,ri) , 
S;{k,v)=Ss{-k,r^). 



(4.41) 
(4.42) 
(4.43) 
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Then, from the definition of the vector source term, we obtain the Fourier amplitudes 
S^{k,r]) = -^^ j £~k ■2ke,y6P^Uk,v)Spiik - ~k,r]) , (4.44) 
S2{k,7]) = [ £k ■ 2kei~k'6PnUk,v)Spi{k ~ ~k,v) , (4-45) 



(27r)3/2 

Ssik,r])=0. (4.46) 

This last equation tells us that the projection of Eq. (4.38) onto the basis vector ki gives 
zero (since contracting Eq. (4.38) with fcj contracts two copies of k with the permutation 
symbol, automatically giving zero). The complex conjugates are 

Slik,r]) = j d^k2ke,~k'6P^Uk,v)Spi-ik + ~k),v) , (4.47) 

S;{k,r]) = j d''k2ke,k'6P^Uk,v)Spi-ik + k),r]) . (4.48) 
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Having now obtained the required amplitudes and their conjugates, we can compute 
the correlator of the source terms for the A = 1 mode in Eq. (4.37).^ Assuming that 
the fluctuations 6p and 5Pnad are Gaussian, we can put the directional dependence into 
Gaussian random variables E{k), which obey the relationships 

(E*(fci)^(fc2)) = 5\k, - fcs) , (^(fci)^(fc2)) = 5%k, + k2) , (4.49) 
and write, for example,^ 

Sp{k,r]) = Sp{k,7])E{k). (4.50) 

The correlator then becomes 

{S*ik,,fj,)S{k,,f]2)) = -^ j d%2hkue\SPr,Uh,Vi)Sp{\k^ + k^\,f],) 

x/d3^,2a,.e-5P..,(^,,.^,)^ 

X (^E{-k^)Ei-k^ - k^)E{k2)E{k2 - h)) , (4.51) 

where we have introduced the notation /j = f{fji)- Wick's theorem (see, e.g., Ref. [16]) 
then allows us to express the correlator in the above in terms of delta functions as 

(^(fci)E(-fci - k,)E{k2)E{k2 - k)) = (^(fci)i(-fci - fci)) (^E{k2)E{k2 - h)) 

+ (E{ki)E{k2)) (^(-fci - ki)E{k2 - h)) 
+ {E{kr)E{k2 - fcs)) (^(-fci - k^)E{k2)) 
= 6%ki + fc2)5'(-fci - fci + fc2 - k) + 5'(fci + k2- k2)6'\-ki - fci + fca) , (4.52) 



■^Note that we need only consider one orthogonal component of the source vector, since we can make 
an appropriate choice such that its component in one direction is zero. 

^In making this choice, we are assuming that Sp and (5Pnad are completely correlated variables. This 
is perhaps not the most physically motivated assumption, since one might expect some level of 
decorrelation between the two variables. However, this assumption will likely give the largest 
signal (the partially decorrelated case will, in its simplest form, require a new parameter less than 
one, which characterises how decorrelated the two variables are - this parameter can be determined 
from the specific model for the production of the non-adiabatic pressure perturbation) , so is suitable 
as a first approximation. We leave the case where the two variables are decorrelated for future work. 
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which gives 

(5*(fci,?7i)5(fc2,r/2)) 

X k2eikl6Pr,^d{k2,V2)Sp{\k2 - fc2|,r?2)| (4.53) 

X + k2)5%k2 -ki-ki- k2) + 5\ki -k2 + k2)5\k2 - fci - fci)} . 

(4.54) 

By integrating over the delta functions, and collecting terms, we arrive at 



{S\k,,fl{)S{k2,fl2)) 



5\k2-k{)k^ / d^^e,k'f5P^^A{k,f|^)5p^{\k + klf|{] 



27r3 

5^nad(|fc + k\,V2)5pi{k,f]2) - 5Pnad (^, ^2)5pi ( | ^ + ^[,^2) 



(4.55) 



from which we can read off the power spectrum: 



47r4 




fihVi V2 df]idT]2 / d k{eik') 6Pnad{k,r]i)6pi{\k + k\,f]i] 



SPns.d{\k + k\,r]2)Spi{k,r]2) - SPni,d{k,r]2)6pi{\k + k\,r]2) . (4.56) 



Now, as a first approximation for the source term, we can expand Eq. (3.42) to lowest 
order in krj to give 



5p{k,r]) = Ak^T]-^ 
SPn,d{k,v)=Dk''V~' 



(4.57) 
(4.58) 



where A and D are yet unspecified amplitudes, and a and (3 undetermined powers. 
Using these approximations gives 



In 



J d^kieik^fk^lk + kf(\k + - + kf 



(4.59) 



where we have performed the temporal integral by noting that as mentioned above, 
a (X 1] during radiation domination, and thus po oc ?7~^. To perform the fc-space 
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integral, we first move to spherical coordinates oriented with the axis in the direction 
of k. Then, denoting the angle between k and A; as 9, the integral can be transformed 

as 

rkc _ _ PIT 

d^~k-^ 271 / Pdk / sin^ci^, (4.60) 
Jo Jo 

where the prefactor comes from the fact that the integrand has no dependence on the 
azimuthal angle, and kc denotes a cut-off on small scales. This cut-off is chosen to be 
smaller than the typical separation of galaxies, and therefore much smaller than the 
continuum limit, solely for the purpose of studying this toy model. Noting that, in this 
coordinate system /cjC* = k sin 6 the integral in Eq. (4.59) becomes 



I{k) = 2TT / k'^+'' sin e de dk sin^ ek^ 1 + {k/kf + 2{k/k) cose 
Jo Jo 

X (k^'ii + {k/kf + 2{k/k) cos ey'^k^ - k'^k^ii + {k/kf + 2{k/k) cos ef/^ ' 

(4.61) 

Finally, we change variables again to dimensionless u and v defined as [11] (or similarly 
[27]) 

V = p u= yjl + {k/kf + 2{k/k) cosO = Vl+v^ + 2v cos 9 , (4.62) 
for which the integral (4.61) becomes 



0/2 



/(A;) = A;^("+«+s /" ' r uduv^dvu^v''(l-^{u'~v'-lf] 
Jo J\v-i\ ^ 4f^ / 



(4.63) 



4.4.2. Evaluating the Vorticity Power Spectrum 



In order to perform the integral Eq. (4.63) derived above, we need to specify the ex- 
ponents for the power spectra of the energy density and the non-adiabatic pressure a 
and p. The energy density perturbation on slices of uniform curvature can be related 
to the curvature perturbation on uniform density hypersurfaces, (, during radiation 
domination through [112] 

^/^ = -^C = 4poC, (4.64) 
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and hence the initial power spectra can be related as (Vini^Pini) = ISpo^ini (CiniCmi), and 
we get the power spectrum of the initial density perturbation 

5 Pirn OC Cini OC f — j , (4.65) 

where fco is the Wmap pivot scale and ng the spectral index of the primordial curvature 
perturbation [TV]. This allows us to relate our ansatz for the density perturbation, to 
the WMAP-data which gives 

5P = Spir^i (^](^] ' = AniPOini ( ^] '^"^'^ (-]' ■ (4.66) 



7/1 I 1111/" Uiill I 7 I I 

koj XVoJ \ko/ \Vo 

Prom this, we can read off that /9 = ^(ng + 1) ~ 1 and the amplitude A = AjniPoini- We 
have some freedom in choosing a, however would expect the non-adiabatic pressure to 
have a blue spectrum, though the calculation demands a ^ /3. Using the notation of 
Ref. [77] we get 

A' = plMkor = pLik^'Ai , = plMko? = pI^^o'^% (4.67) 

where we also have the ratio^ 

A| 1 - a(fco) ^ ' 

and therefore, 

{ADf = -^^-^^'nPLik^'' ■ (4.69) 

We can then substitute in numerical values for and a(/co) from Ref. [77] later on. 
Then, making the choice a = 2, the input power spectra are 

SpiK n) = (^) (^) . ^^'-C^. ") = {£) " . (4-0) 

for which the integral Eq. (4.63) becomes 

kc/k I'V+l 1 



I{k) = k'' I ^ u'duv'^dv^l- ^{u' -v' -ly) u'v-v'u . (4.71) 



^The parameter a(fco) is introduced in Refs. [22, 77] in order to quantify the ratio of to A|j. 
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We can then integrate this analytically to give 



135 ' 245 ' 1575 ' 



(4.72) 



which clearly depends upon the small scale cut-off, as expected.'' For illustrative pur- 
poses, we choose kc = lOMpc"^ and plot the solution I{k). Fig. 4.1 shows that the 
amplitude of the integral grows as the wavenumber increases. Fig. 4.2 shows a turn 
around and a decrease in power at some wavenumber (in fact, for a non-specific cutoff, 
this point is at 3.7375/cc). However, we note that this value is greater than our cutoff, 
and therefore not physical. As long as we consider values of k less than the cutoff, our 
approximation will still be valid. 

Then, using the above, and noting that the input to the temporal integrals are 



a (X 



Po = POii 



(4.73) 



we obtain the power spectrum for the vorticity, for general kc as 

2 



:kn 



1 



256 '^47rVLi 



In^ 



V 



a{kt 



0) 



1 - a{ko) 



. -12 A 8 4 , 

'^0 '-^nPoim'^- 



16 kf: ( k 



135 fcg V^o 



+ 



12 kl ( k 



245 A;^ V^o 



1575 Vk^ 



11 



(4.74) 



Substituting in values of parameters from Ref. [77], as presented in Table 4.1 taking a 
conservative estimate for a(fco), being 10% of the upper bound as reported by Wmap7, 
we obtain the vorticity power spectrum 



Vo 



0.87 X lO^^^A;? 



ko 



Mpc" + 3.73 X lO-'^^fc, 



18^7 



k 



7.71 X lO^^^A;,^ ( ^ ) Mpc^ 



11 



Mpc^ 



(4.75) 



^It should be noted that if one were to reduce the assumption of 100% correlation between the energy 
density and entropy perturbation, this could soften this dependence. However, this is left for future 
investigation. 
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Figure 4.2.: Plot of I{k) for the illustrative choice of = lOMpc ^; wide range of k 
including k > kc- 
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Parameter WMAP7 value 



ko 0.002 Mpc"^ 

A^{ko) 2.38 X 10"^ 
alko) 0.13 (95% CL) 



Table 4.1.: Parameter values from the WMAP seven year data [77]. 



and for our above choice of kc = lOMpc ^, 
r4k,r])=r]Hn'' 



0.87 X 10 



-3 



kV 



+ 3.73 X 10 



-11 



-7.71 X 10 



-20 



kn 



11 



Mpc^ 



(4.76) 



This shows that, under our approximations, the vorticity spectrum has a non-negligible 
amplitude, with a huge amplification of power on small scales. We plot this, for illus- 
trative purposes, in Figs. 4.3 and 4.4, where we have ignored the time dependence and 
focused only on the scale dependence of the spectrum. We should emphasise that we 
are studying the generation of vorticity in the wavenumber region k^ < k < k^ and do 
not expect the dynamics to be dominated by an "inverse cascade" (i.e. a feedback of 
power from smaller scales to larger scales). Therefore, the physics around the cut-off 
wavenumber cannot influence the vorticity generation. To study this phenomenon in 
detail, a much more detailed calculation including backreaction effects would have to 
be performed, beyond the scope of this thesis. 



4.5. Discussion 



In this Chapter we have studied the generation of vorticity in the early universe, showing 
that second order in cosmological perturbation theory vorticity is sourced by first order 
scalar and vector perturbations for a perfect fluid. This is an extension of Crocco's 
theorem to an expanding, dynamical background, namely, a FRW universe. Whereas 
previous works assumed barotropic flows, allowing for entropy gives a qualitatively 
novel result. This implies that the description of the cosmic fluid as a potential flow, 
which works exceptionally well at first order in the perturbations, will break down at 
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Figure 4.3.: Plot of Vu>(^), i-e. the scale dependence of the vorticity power spectrum. 
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second order for non-barotropic flows. Similarly, in barotropic flow Kelvin's theorem 
guarantees conservation of vorticity. This is no longer true if the flow is non-barotropic. 

Having derived the qualitative result, we then obtained the first realistic calculation 
of the amount of vorticity generated at second order. As an input spectrum for the 
linear energy density, we used the solution obtained in Section 3.1.2 approximated for 
small krj and normalised to Wmap7. Then, making the ansatz that the non-adiabatic 
pressure perturbation has a bluer spectrum than that of the energy density in order to 
keep the non-adiabatic pressure sub-dominant on all scales, we obtained an analytical 
result for the vorticity. Our results show that the vorticity power spectrum has a non- 
negligible magnitude which depends on the cutoff, kc, and the chosen parameters. As 
this is a second order effect, the magnitude is somewhat surprising. We have also shown 
that the result has a dependence of the wavenumber to the power of at least seven for 
the choice a = 2, where a is the exponent of the wavenumber for the non-adiabatic 
pressure input spectrum. Therefore the amplification due to the large power of k is 
huge, rendering the vorticity not only possibly observable, but also important for the 
general understanding of the physical processes taking place in the early universe. 

The consequences of this significant power are not immediately clear as the model 
under consideration is only a toy model. Although vorticity is not generated in standard 
cosmology (at linear order), the vorticity generated at second order will not invalidate 
the standard predictions, as it is, on large scales, very small and can only be of significant 
size on small scales. However, any possible observational consequences will depend on 
the wavenumber at which the power spectrum peaks (to be determined by an actual 
model). 

One prospect for observing early universe vorticity is in the B-mode polarisation of the 
CMB. Both vector and tensor perturbations produce B-mode polarisation, but at linear 
order such vector modes decay with the expansion of the universe. However, vector 
modes produced by gradients in energy density and entropy perturbations, such as those 
discussed in this chapter, will source B-mode polarisation at second order. Furthermore, 
it has recently been noted that vector perturbations in fact generate a stronger B-mode 
polarisation than tensor modes with the same amplitude [55]. Therefore, it is feasible 
for vorticity to be observed by future surveys such as the space-based CMBPoL [21], 
which is currently in the planning stage, or the ground based experiment POLARBEAR 
[53], which is due to start making observations in 2011. 

Finally, a non-zero vorticity at second order in perturbation theory has important 
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consequences for the generation of magnetic fields, as it has been long known that 
vorticity and magnetic fields are closely related (see Refs. [24, 61]). Previous works 
either used momentum exchange between multiple fluids to generate vorticity, as in 
Refs. [58, 72, 74, 103, 115, 139, 146], or used intermediate steps to first generate vorticity 
for example by using shock fronts as in Ref . [ I "^^ •] . However, we do not require such 
additional steps. Therefore, an important extension to the work presented in this 
chapter is to consider the magnetic fields generated by our mechanism which could 
be an important step in answering the question regarding the origin of the primordial 
magnetic field. We will discuss more future prospects in the concluding chapter of this 
thesis. 



5. Third Order Perturbations 



In the preceding chapters we have developed cosmological perturbation theory at both 
hnear and second order. However, one does not need to stop there: it is worthwhile 
and feasible to consider perturbation theory even beyond second order. In this chapter 
we explore aspects of cosmological perturbation theory at third order. 

We have already seen that extending perturbation theory to second order reveals 
new phenomena which arise due to quadratic source terms in general, and the coupling 
between different types of perturbation which is not present at linear order. At third 
order a new coupling occurs in the energy conservation equation, namely the coupling 
of scalar perturbations to tensor perturbations. This will allow for the calculation of yet 
another different observational signature, highlighting another aspect of the underlying 
full theory. 

There has already been some work on third order theory. For example, Refs. [69, 71] 
considered third order perturbations of pressureless irrotational fluids as "pure" general 
relativistic correction terms to second order quantities. The calculations focused on the 
temporal comoving gauge, allowing the authors to consider only second order geomet- 
ric and energy-momentum components, and neglected vector perturbations. Ref. [85] 
includes a study of third order perturbations with application to the trispectrum in 
the two-field ekpyrotic scenario in the large scale limit. There has also been reference 
in the literature of the need to extend perturbation theory beyond second order. For 
example, in Ref. [ • '] UV divergences in the Raychaudhuri equation are found when 
considering backreaction from averaging perturbations to second order. The authors 
state that these divergences may be removed by extending perturbation theory to third, 
or higher, orders. 

In this chapter, we develop the essential tools for third order perturbation theory, 
such as the gauge transformation rules for different types of perturbation, and construct 
gauge invariant quantities at third order. We consider perfect fluids with non-zero 
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pressure, including all types of perturbation, namely, scalar, vector and tensor pertur- 
bations. In particular allowing for vector perturbations is crucial for realistic higher 
order studies, since vorticity is generated at second order in all models employing non- 
barotropic fluids as shown in Chapter 4. Hence studying irrotational fluids at higher 
order will only give partial insight into the underlying physics. We present the energy 
and momentum conservation equations for such a fluid, and also give the components of 
the perturbed Einstein tensor, up to third order. All equations are given without fixing 
a gauge. We also give the Klein-Gordon equation for a scalar field minimally coupled 
to gravity at third order in cosmological perturbation theory. This work is published 
in Ref. [3(J]. 



5.1. Definitions 

As in Chapter 2 we take the perturbed metric tensor with covariant components 

5foo = -a^(l + 20) , goi = a^Bi, gij = {6ij + 2Cij) , (5.1) 

In order to obtain the contravariant metric components we impose the constraint 
Eq. (2.14), g^vQ^'^ = 5^'^, up to third order. This gives 

(7°° = (l - 20 + 402 - 803 - BkB'' + ic^BkB^ + 2B'B^Ci^ , (5.2) 

g^i = 1-[b'- 2(t>B' - 2BkC^' + A(fB' + ABkC^'cp + AC^^C^'B^ - B'^B^B'^ , (5.3) 
gij = 1 f^ij _ 2C'^ + AC'^Ck^ - B'B^ + 2(j)B'B^ - ^C'^C^^Cki 

+ 2B'C''^Bk + 2BkB^C''^^ . (5.4) 

Note that in this chapter we do not explicitly split terms up into first, second and third 
order parts unless where necessary, since doing so would dramatically increase the size 
of the equations presented. For example, expanding the — component fully order 
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by order gives 

g^^ = -- (l- 201 - 02 - ^03 + 02 + 80102 - 80? - fii,fii*= - B^kBi' 



+ A<P,B,^B,k + 2B:B,^Cu, , (5.5) 



an 3 

40i5i^5i, + 2B,'B,^Cuj J , 

which when compared to Eq. (5.2) illustrates the increase in number of terms, and thus 
why we refrain from splitting perturbations up. 

Furthermore, to third order in perturbation theory, the fluid four velocity, defined in 
Eq. (2.21) and satisfying the constraint 

M% = -1, 

has components 

(5.6) 

+ ^Vkv'' + VkB^ + CkjV^v^ - 2^v''Bk - (pv'^Vk^ , (5.7) 
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a 


u' = 
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Ui = 


a (^Vi + 5j - 


Uo = 


-a ( 1 + - 



-02 - -^^ 

2^ 2 



05, + 2akv'' + ^S,02 + ^B.v'^Vk + B.v^Bk 



1 



(5.8) 
(5.9) 



5.2. Gauge Transformations 

We firstly need to extend the gauge transformations derived earlier to third order in 
cosmological perturbation theory. Expanding the exponential map, (2.50) 

T = e^«T, (5.10) 

to third order gives 

exp(i:5) = l + t£^, + ]^^£l + ]^^£i:, + ^tH^, + \(^''£l. + -/£ii£^2 + \(^''£i2£ii. + ■■■ 

(5.11) 
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Splitting the tensor T up to third order and collecting terms of like order in e we find 
that tensorial quantities transform at third order as 

+ 3{£l + £^,)6Ti + 3£^,6T2. (5.12) 

Then, by expanding the coordinate transformation, Eq. (2.56) to third order we 
obtain the relationship 

(5.13) 



1 3 



relating the coordinates of the points p and q. 

In the rest of this section, we will derive the gauge transformations at third order in 
an analogous way to those derived for first and second order in section 2.3. 

5.2.1. Four Scalars 

At third order we also split the generating vector ^3 into a scalar temporal and scalar 
and vector spatial part, as 

^3^= («2,/32; + 730 ' (5-14) 

where the vector part is again divergence-free {dk'^^^ = 0). We then find from Eqs. (5.12) 
that the energy density transforms as 

+ 3 + £52) 5pi + 3i'5^(5p2 , (5.15) 
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which gives 

(^Ps = (^Ps + Potts + Po'tti + 3poaiai,A^i + Po {aiMB^i + ai,x^i, /?) 
+ Spo'aias + Po- (tts.A-^i + tti^A^s) + 3 {Spi,xp^i + Spi^x^l ^) 
+ 3(5^1,^^2 + 3(5p2,Aei • (5.16) 

Similar to the second order case, we need to specify the time shcings (at all orders), 
and also the spatial gauge or threading at first and second order, in order to render the 
third order density perturbation gauge- invariant. 

5.2.2. The Metric Tensor 

We now give the transformation behaviour of the metric tensor at third order. The 
starting point is again the Lie derivative, which for a covariant tensor is given by 
Eq. (2.68). 

As above in the case of the transformation behaviour of a four scalar at third order, 
the change under a gauge transformation of a two-tensor can be found applying the 
same methods as at second order. We therefore find that the metric tensor transforms 
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at third order, from Eqs. (5.12) and (2.68), as 



,(0).A 



o„(0) , (0) /"tA t/3 , fA , ^ (0) /.A f/3 , tA t/J 

+ (5'1a,/3'^2,i/ + 9x^,i3^2,fi + 9ltX^2,l3u + 9x,!^2,l3fij ^1 + (^fi'/jj,A'^2 + 9xi!^2,l3^ ^1, 
+ (fi'^J,A^2 + 9xi!^2,l3^ + (^fi'^A,/3'Cl,iy + 9xi!,l3^1,fj. + 9\2^1^pv + 9xi!^l,l3fij C 

+ (5'^J,A^1 + fi'li^l,/?) (s'^i.A^l + 9xi!^l,P^ ^2, 

I J ^(0) faf/3(^A , 0^(0) f AfOf/? , ^(0) , ^(0) ^A f/3 

"T ^ i//x!/,a^ASl SI SI "T -jy^i/.A/gSlSl Sl,a "T y^iy,ASl,a/3Sl SI "I" i//i;/,ASl,/3Sl,aSl 

"T ^9fj.p^X Sl,aSl,i/ "T Si Si S1,q;/ y "T 'Jy/ia.ASl Sl,^Sl,i/ "r '-'i//ia,/3ASl SI Sl,j/ 



, (0) 



^OLf^l^/^X I t"A I /^a/^f^ ^X I r\/:0/^X ^OL _i_ ^(^ ^X /-a 

SI SI Sl.a^j/ "T Sl,aSl,/3Sl,jy "r SI Sl,aSl,/3jy "r ^Sl Sl,a/3Sl,i/ "t" SI Sl,oSl,/3v 



"T ^9ul3,X I SI Sl,aSl,^ "I" SI SI Sl,Q/^ I "T ■^9^a,X^l Sl,/3Sl,/^ "I" '^9^a,l3X^l SI Sl,^ 



(0) 
^.A 



tatfi/^X I t-A _i_ catfi /^X , n/:l3^X ta , cf^cX ca 

SI SI Sl,Q^^ "T Sl,aSl,/3Sl,^ "T Si Sl,aSl,^/x "^Sl Sl,a^Sl,^ "r Si Sl,cjSl,/3^ 



tA / f a f /3 _|_ fCK f/3 \ I f/3 /' fA i fo: C-^ ^ 
Sl,/3 I Sl,^Sl,j/ "I" Sl,i/Sl,^ / ~'~ "^1 l^l.At^l,/?'^ ~'~ 



(5.17) 



However, in this case it becomes even more obvious than in Section 5.1 that the ex- 
pressions at third order are of not inconsiderable size. This will also be clear from the 
Einstein tensor components and the evolution equations given below in Section 5.4. 
Now, following along the same lines as at second order, Eq. (5.17) gives the transfor- 
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mation for the spatial part of the metric at third order, 



(5.18) 



where X^ij contains terms cubic in the first order perturbations. Extracting the curva- 
ture perturbation gives 



3 ,ij 



(5.19) 



This expression is general, including scalar, vector, and tensor perturbations and is 
valid on all scales. However, we shall detail here only the expression valid for scalar 
perturbations and large scales and find that X^ij takes then the simple form 



X: 



3ij 



+ 2na2ipi + naiip2 + 2\^—+ W j afipi 



+ 



"H (aiCCg + a[a2) + 2 ( h ) aia2 



+ 3n- 



(5.20) 



Hence we finally get for the transformation of ■03 



--03 = --03 + + ( — + Sn— ]al + 3l—+n'] aia[ 

a a J \ a 



-3 



(^2ip'i + 2 "1^2 + ii^'i + 2^0i) + ctl « + mij'i) + 2Ha2ipi 



a 



+ 'Haiij2 + 2[— + W] aiipi 



(5.21) 



5.3. Gauge Invariant Variables 



In the previous section we have described how perturbations transform under a gauge 
shift. We can now use these results to construct gauge-invariant quantities, in particular 
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the curvature perturbation on uniform density hyper surf aces, (. In this section, as be- 
fore, we consider only scalar perturbations, and restrict ourselves to the large scale limit. 

We first define hypersurfaces in different gauges as in Section 2.4. From Eq. (5.16) 
we find the time slicing defining uniform density hypersurfaces at third order in the 
large scale limit as 



Sp3 ^ _J_ 



3 (5pi5p2 + Vl^P2) 



6p'l6p ' 



1 ,2Spi II I 

Po Po 



5pi 
Po 



(5.22) 



Similarly, the temporal gauge transformation on uniform curvature hypersurfaces is 
defined by evaluating (5.21) and gives, at third order. 



ttSflat 



•H 2?/2 



n 



37, 



8ipf 

~7r 



(5.23) 



We can now combine the results found so far to get gauge invariant quantities, and 
as before choose the curvature perturbation on uniform density slices as well as the 
density perturbation on uniform curvature hypersurfaces as examples. 

One gauge invariant matter quantity of interest is the perturbation to the energy den- 
sity on uniform curvature hypersurfaces. This is obtained by substituting the temporal 
gauge transformation components in the uniform density gauge into the appropriate 
transformation equation, Eq. (5.16). This gives, at third order. 



5p: 



'3flat 



+ 3 



Sp3 + 



/ ^3 



1^(2^2^; + V^^V^i) + 3^ [2(p>i + ^IPo) + < 



Po + 2po^ 



+ 



/// 

Po 



3Po 



Po- 



m 



Po 



. (5.24) 



The curvature perturbation on uniform density hypersurfaces, (, as introduced in 
Eq. (2.126), is defined as 

-C = i^p- (5.25) 
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This is obtained by substituting the temporal gauge transformation components in 
the uniform curvature gauge into the appropriate transformation equation, Eq. (5.21). 
Evaluating this on spatially flat hypersurfaces then gives, at third order 



Po Po 



-;^Sp26pi 
Po 



-yr-Opi dpi j^dpi dp 

Po Po 
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Po 
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a a 



3^(-+W 
Po V a 



(5.26) 



There are different definitions of the curvature perturbation present in the literature, 
depending on different decompositions of the spatial part of the metric tensor. A differ- 
ent definition to the one above, as discussed in e.g. Ref. [H-], was used by Maldacena 
in Ref. [10 i] to calculate the non-gaussianity from single field inflation, and was intro- 
duced by Salopek et al. in Refs. [130, 131]. They define the local scale factor a = e°, 
then 



,2o 



Comparing to the expansion from perturbation theory 



(5.27) 



(5.28) 



one can obtain the relationship 



C - CSB + CsB + gCsB • 



Splitting this up order by order gives, at second order 



(5.29) 



C2SB — C2 — 2(Ci 



(5.30) 



and, at third order, 



CasB = Cs - 6C2C1 + 8(Ci 



(5.31) 



Note that it is this definition, (5.31), of the curvature perturbation which occurs in 
Ref. [85] , though with different pre-factors since their perturbative expansion is defined 
differently. It is perhaps worth mentioning that Csb, the variable first introduced by 
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Salopek and Bond and then employed for non-gaussianity calculations by Maldacena, 
is extremely Gaussian after slow-roll inflation, as opposed to other ( variables which 
exhibit non-gaussianity, as can be seen e.g. from Eq. (5.30). 

5.4. Governing Equations 

Having constructed gauge invariant quantities up to third order in the previous section, 
we now turn to the evolution and the field equations. The equations presented in this 
section in full generality are new; Refs. [69, 71] previously considered some governing 
equations at third order in perturbation theory, however they focussed on pressureless, 
irrotational fluids. 

5.4.1. Fluid Conservation Equation 

In this section, we give the energy momentum conservation equations for a fluid with 
non-zero pressure and in the presence of scalar, vector and tensor perturbations. The 
latter generalisation is important since at orders above linear order, all types of pertur- 
bation are coupled. 

As at linear and second order, presented in the previous chapters, energy momentum 
conservation 

V^T^ = 0, (5.32) 

gives us evolution equations. Substituting the definition for the energy momentum 
tensor, Eq. (2.18) expanded to third order, into Eq. (5.32) gives energy conservation 
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(the 0-component) 



V + m{5p + 5P) + (po + Pq){C% + v,^) + {5p + 5P){C% + t;^;) + ((5p + 5P)^iv' 
+ (po + Po) [(5' + 2i;^)(i;: + B[) + - 2C;.C*^' + v\C^ + 20,,) + 4Ht;'(5, + 2i;. 
+ (5p + 5P) \{B' + 2i;')(t;: + B[) - 1C[p' + v\C^ + 20,, + 4H(t;i + 5,)) + t;^^. 
+ (5p' + 6P'y{Bi + t;,) + ((5p + 6P)^i(Pv' + (po + Po)i2C'^ViV, - B^v'^ + v'v^) 



- (po + Po) 2a%C,,,kV^ - 2v[v, + B,Br - ^C^^C]^) + -v\{<^^ - v'v,) 



+ [v' + B'){2B[(t) - C^Vi) - vHB'Bij + v'vij + (f){vj + 2v'j + SUvj - 2(f) j + C\j 



- 2(1)' B J + 3C',/ \ + B^{B'C[^ + B,(P' + v'^ 



(5.33) 



and momentum conservation (the i— component) 

(po + Po)ivi + 5i)]' + (po + Po)(0,i + 4:n{v, + Bi)) + 5P,, + [(5p + 5P){vi + 5^) 
+ (5p + 5P)(0,, + 47{(t;i + B,)) - (p'o + P^) [(25^ + t;,)0 - 2Ci,t;^' 
+ (po + Po) [(t^i + Bi){0'^ + t;^',,) - B,{(t)' + 8H0) + i;^(Pi,, - P,,, + t;,,,- + SHC,, 
+ (2Q,t;^y - 0(t;^ + 2P,' + 20,, + AV,u 



-2C,,0) + (-t;, + 4P,)02 



+ (Po + Po) p(Pi^j + Pj^i + 2P,Pj + -WiWj 
+ (po + Po) 0' ( 4P,' + + 2H(8Pi + 3^;^) + 4(/ 



+ (t;, + 2Bi){v'^B^ - C^;.0) + (t;, + B,){B'^B^ - C'^.C^') + [2B,{B'. + v'^) 

+ P,(2P,' + 20,, + v'i) + 2V,{vi + 2Pi)(2P,- + Vj) + 2C^,{C\ + t;^^ - 4-H0) 
+ C\,j(Pi + t;,) + Vi{B'^ + 0,j) + Vj{B[ + 0,,) + w'^(2C,fc,j. - Cjk^i) 

+ (P,,, - P,,, - 2C:,)0 + 2C,fci;^,| + \{v,v,v^)' - 2C,,v^'(P + P,(40'0 - v^^ 



= 0. 
(5.34) 

As emphasised earher, Eq. (5.33) highhghts the couphng between tensor and scalar 
perturbations which occurs only at third order (and higher) in perturbation theory. At 
both linear and second order, no such coupling exists, since the only terms coupling 
the spatial metric perturbation, C,j, to scalar perturbations contain either the trace or 
the divergence of Cij and the tensor perturbation, hij is, by definition, transverse and 
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trace-free. However, at third order, terms like 5pC[jC'''^ occur in the energy conservation 
equation which, on sphtting up order by order and decomposing Cij becomes 6pih[-jhi''K 
It is clear that this term only shows up at third order and beyond. Thus, as mentioned 
earlier, third order is the lowest order at which all the different types of perturbations 
couple to one another in the evolution equations, which will produce another physical 
signature of the full theory. 

Scalars only 

It will be useful to have energy and momentum conservation equations for only scalar 
perturbations. These equations are obtained by making the appropriate substitutions 
Cij — —Sijip + E^ij, Vi — v^i, and Bi — Bi into the above expressions. On doing so, we 
obtain the energy conservation equation 



5p' + m{5p + 5P) + (po + i^o)(V\ + V'^E - 3V^') + {5p + 5P){V''v + V'^E - 3V^') 



+ {5p + 5P)^iv: + (po + Po) [{B' + 2v:){v[, + s;.) + - 2 (^''(3^' - V'^;) 



+ {5p + 5P) [{B: + 2v:) {v[, + B[,) - 2 (^''(SV' - V^^;) - i^V'E' + E^E/^) 
+ v^{2(t),i - 2>ij^i + V'^E + m{v^i + + V''v<p\ + (V + 5P')v^{B^i + v,i) 
+ {5p + 5P),(t>v: + (po + Po){v:v,i(t> - B,iv:<P - G^jv^^v, + E:^v,v,,) 



- (po + Po) [2E:^{v,''E,ijk - 2v[,vj + B,iB[^ - 2E/E,kj) - 2V'£;V,fe^,' 

- 2ij(v^\V^E,k - 3^,k) - 2t;>/ + B[^B/ + 6^^') + ^V'w((/.' - vfvj) 




(5.35) 



and the momentum conservation equation 
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+ [{5p + 5P){v^i + S ,)] ' + (5p + (5P)((/),, + 4?{(^;,, + S ,)) 

- (p[, + i^) [(25,, + t- + 2^t;,i - 2E,,t;/ 

+ (po + Po) [{v,i + 5,i)(V2^; - 3V^' + V^^;') - P,,(0' + m(t>) 

+ - mitl^Sij - - 2(V'T;,i + - + 2P;, + 20,, + 47^^;,,)] 

+ (Po + ^o) [^,KB,^,j + + 2P,,P,,- + ]^v,^,^ + 20(V^5,,- - E,,)) 

+ + 4P, + (po + Po) + 2P,,) ^' - (3V^' - V'£;')0) 

+ (^4P;, + ^-v[, + 2?{(8P , + 3^;,,) + 40,,^ + {v,^ + P ,) (p;,.P - 3^^^ 

+ ^'v'p + ^v^'p' - p/''p;,fe) + t;,^{2P,,(p;^. + 1;;^.) + p,,(2p^ + 20,i + v[;) 

+ 2?^(^,, + 2P,,)(2P,,. + v,^) - (3^^,,- - V2p,-)(^,i + ^,0 

- 2(V'5i,- - E^ij){V\ + V'P' - 3^^' - 47^0) + ^;,,(P,'^. + 0,,) 

+ ^^(P^ + (/),i) - 2V'j'y,i + Pifej'y,'' + i^,iv,j + 2(V''5y - Py))0 - 2#,ij- 

+ 2Pife^;,^} + \{v,iV,jV,')' + 2^^;;,0 - 2E,ijV,''<i) + Pi(40> - V'l;^)] = . (5.36) 

Considering the large scale limit, in which spatial gradients vanish, the energy con- 
servation equation becomes 



V + ?>U{5p + 5P) - 3V''(po + Po) - 3V''((^P + 5P) 

- 6VV'(po + Po) - {5p + 5P) + 12V''V'(Po + Po 



(5.37) 



Splitting up perturbations order by order, this becomes 



5p^ + m{5p^ + 5P3) - 3V'^(po + Po) - 9V'^(5pi + 5Pi) - 9V';((5p2 + 5P2) 

- 18(po + Po)(V'2V'i + V'iV'2) + 72V'iVi(Po + Po) = . (5.38) 
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In the uniform curvature gauge, where ip = 0, this is 

Vsflat + 3?/(5p3flat + ^Aflat) = , (5.39) 

and in the uniform density gauge, where 6p = 0, 

SnSPssp + S^ipo + Po) + X'^SPisp + 9C[5P2Sp 

- 18(po + Po)(C2C; + C1C2) - 72Ci'C;(Po + ^o) = O , (5.40) 

with ( as defined above. This can be recast in the more famihar form by introducing 
the (gauge invariant) non-adiabatic pressure perturbation. At hnear order the pressure 
perturbation can be expanded as, from Eq. (2.149), 

dP dP 

= 1^^^^ + l^^Pi = -^^nadi + c^Spi . (5.41) 
ob op 

This can be extended to second order [37] and higher by simply not truncating the 
Taylor series: 

5Pnad2 = SP2 - Cl5p2 ' ^6p,' , (5.42) 

op 

Thus, in the uniform density gauge, the pressure perturbation is equal to the non- 
adiabatic pressure perturbation at all orders. Then, Eq. (5.40) becomes 

C3 + —^5P^.dS = 6(C2C( + C1C2) + 24Ci'C; ^(Ca^i^nadl + ClSPn.d2) ■ (5.44) 

Po + ^0 Po + ^0 

In the case of a vanishing non-adiabatic pressure perturbation, ([ and Q are zero and 
hence we see that ^3 is also conserved, on large scales. This was also found in Ref. [85], 
and previously in Ref. [12]. 

5.4.2. Klein-Gordon Equation 

The energy momentum tensor for a canonical scalar field minimally coupled to gravity 
is easily obtained by treating the scalar field as a perfect fluid with energy-momentum 
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tensor (c.f. Chapter 2) 

T^u = g''^^,x^,u - 5^ (^^g'"'^,a^,f, + U{v)^ , (5.45) 
where the scalar field (f is split to third order as 

^iv, x') = Mv) + ^MV, x') + ^S^2{v, x') + ^S^siv, x') , (5.46) 
and the potential U similarly as 

U{ip) = Uo + 6Ui + ^6U2 + , (5.47) 

where we define 

6U3 = U^^^^5ip\ + 2U^^^5ifi5ip2 + U^^6(p3 , (5.48) 
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and making use of the shorthand notation U^^p = Then, Eq. (5.32) gives the 
Klein-Gordon equation 



6(5 tpo 



ifi[^i2(t>-4:(t>^+B'Bi) 



+8<;i¥'o(W-2'H</>-2(/.0+B'('^J,B^+2Vi,i+4WBO+M)|^'.'+<^'-<^'i(l-2<^) 



+S<pi,iiB''+2HB'+2C'^ 



-3S<fi2,i(B"+C^ j;+(f>,'-2C'\j+2HB'] 



B"+C^j^^-2C'\j+<p,^-\ l^+2'HB^-2BjC'^'+B^C^r+B^ Bj,^-B^ B\j 



-B^Bi ,j-2C^^B'.-2C'W''k,j+'iC^^Cjk,''-2C'^<l>,j-2B^'(f>-B^(l)'-2(f),^<l>-mBjC^^ 



-mB^cj>+AC^Wj\k-'2C^^Ck\j -12B\Sip[. + ^5ip'2-)+2A5ip[.{B',p+C'^Bj) 
+6<5(p2,ijC'J-6<5(pi,ij(4C'=JCfc*-B*BJV24W^;,(/)(l-20+4(/i2_^/^_30/^) 

-6v3[|' ^2<t>{l-2^+44>'^)+B^Bi-4B^B,(t)-2B^BiCij^ -6C'^j¥'^(Bi0-2_Bi) 

+6C*;i^{,(l-2<^+4</-2-BJBj.)_6C^j,i</pJ,(B"-2C*'=Bfc-2BV)-6C'»JVo(2<^y 



-4C,j</.-4Cfc,C^)-12^;,C»^' 
+2C,kB'' j+2Bj^,<j>-mB,Bj 



2BiCjk,'^+^B''Ckii — Bj^i—BiB'-+Bi(l>i—BkCj- 



B\i+B^B[-B^B,B\j-B'-BjB\i-iB^B[. 



+2'HB^B,-2B^Bi<j>'-8-HB^B,^+4:B^<j>^,-2B^ i<i>+4B^ i(t>^-B^<j>^, 



+6l/^a2=0. 



(5.49) 



One can again see the couphng between first order tensor and scalar perturbations. For 
example, the 6ipi^iC'^^(f)j contains a term that looks like 5ipi^ihi^ (pi^j, which occurs only 
at third order and beyond. 

Again, we refrain from splitting up the perturbations order by order for ease of 
presentation. Once split up, one can then replace the metric perturbations by using the 
appropriate order field equations. We present the Einstein tensor at third order in the 
next section. Note also that Eq. (5.49) implicitly contains the Klein-Gordon equations 
at first and second order. We refer the reader to, for example, Ref. [iO(j], for a detailed 
exposition of the second order Klein-Gordon equation. 
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5.4.3. Einstein Tensor 

The Einstein tensor, which describes the geometry of the universe, is defined (as shown 
in Section 1.1) as 

G^ = i?^-^5^i?, (5.50) 

where R^^y is the Ricci curvature tensor and R is the Ricci scalar. Here, we give the 
components of the Einstein tensor up to third order: 



3,^-C[j.' + 1{^'^B,-Bj^,^)+2'H{C^ J,^-2C,J,^-4>,^) 



k,ij ■ 



IBj4B\^ +B\')-3H'^{4<P^-B,B')-^B\,Bj^^ -4:nB\,<t,+G°o , 



(5.51) 



a2GO,=C'='^^_,-C,'^,/=-|(Bfc,,'=-V2B,)-2W0,,+8-H0,,</.+C^(2C'=^fc-C*fc,^^ 

+2C''^ I^C^^ — Cj,._^ + |(-Bfc,ij— -Bi,fej)j +B^ {Ckj,i' —C^ k,ij+Cik.k-' —'^^Cij—2'HBj,i^ 
-l(^B,,^+Bj,,^ 4^,' HB^,J~BJ^,){^C'=^:,' -C^'' ,k)-C,k.j{B'' -B^ ,'')+B\,<f,,i 
+</.[B,,^^-V2B,+2(c^ .'-C^'j,,)]-C'=j'Cfcj,,+GOi , (5.52) 
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-\'^C"j+CjiJ''~n(B\j+Bj 



+'5^■{(«'-^)(l-20)+2W(B^fe-C^'='fc+<^')+B'''fc-C''■^fc^-C'="fe+V2(</,+C^)} 



+{c'=;-<^'-B^fc){c»;-l(B^,+i?,/))+c''='(i?,,fc-2c;^.)+c;^-B^'=+0.>,^ 

+{B'='-2C'=^+CS,'=+0,'=)(C,fc/+C\,j-C',,fc) + |B'(Bfc,/-V2B,+4W</.j-2C'=\.j^ '=) 



h{B^k+K.,)~c'^j+<t'.ik-Cki,j'-C,i,k'+^^Ckj+C'i.,k+H[B,,k+Bk.,-2C',^ 



-\{Bk:B^,j+Bj,^B\k)+'t> 

(402_Bj^Bfc)+2</, 



C"'"fc-B'''fc-V20+2W(C'='fc-20'-Bfc,fc) 



2C'',,fe-2C^,/+V2iJfc-B,,fe'+2-H(B^-</.,fc-2C'fc,,+C^,fc) 

Cf,i—'V^Ckl+2'HCf,i+2Clm,k^^ — C'^rn,kl — '^'HBl,k — Bij. — 4>^l^.l 



C^^'ilc'-B, 



+2B'''iC'i,k-Cki,') 



\ r^k f Jul \ ] o^kl m \ ( r^km I i r^lm k\ /^l f or^k m 1 /^m k \ 

+ fcl-O,!— 2^ ;)+^*-' .k'-'lm, +'-^lm.,k['-^ > ~ 2^ ' I" <-'i,fcM<^ m, — 



+<A'(c*^-B^fe)-|(2B^fcB^/-i?^,fcB^'-3B^fci?^,'■■)+</.,4c'^.''-2C'^^-9i,''■)+G^J}, (5.53) 

where G°o, G^j, are the third order corrections (the latter spht into a diagonal part 
Glj, and an off diagonal part GJ^j) which we give in the appendix as Eqs. (B.l), (B.2), 
(B.3) and (B.4), respectively. Note that, in calculating the third order components given 
above, we have implicitly obtained the full second order Einstein tensor components for 
fully general perturbations (i.e. including all scalar, vector and tensor perturbations). 



5.5. Discussion 

In this chapter we have developed the essential tools for cosmological perturbation the- 
ory at third order. Starting with the definition of the active gauge transformation we 
have extended the work presented in Section 2.3 to third order, and derived gauge in- 
variant variables, namely the curvature perturbation on uniform density hyp ersurf aces, 
Cs, and the density perturbation on uniform curvature hypersurfaces. We also relate 
the curvature perturbation ^3, obtained using the spatial metric split of Ref. [112] to 
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that introduced by Salopek and Bond [130], which is also popular at higher order. 

We have then calculated the energy and momentum conservation equations for a 
general perfect fluid at third order, including all scalar, vector and tensor perturbations. 
The Klein-Gordon equation for a canonical scalar field minimally coupled to gravity is 
also presented. We highlight the coupling in these conservation equations between scalar 
and tensor perturbations which only occurs at third order and above. Finally, we have 
presented the Einstein tensor components to third order. No large scale approximation 
is employed for the tensor components or the conservation equations. All equations are 
given without specifying a particular gauge, and can therefore immediately be rewritten 
in whatever choice of gauge is desired. However, as examples to illustrate possible gauge 
choices, we give the energy conservation equation on large scales (and only allowing for 
scalar perturbations) in the fiat and the uniform density gauge. This gives an evolution 
equation for the curvature perturbation ^3, Eq. (5.44). As might be expected from 
fully non-linear calculations [KHl] and second order perturbative calculations [110], the 
curvature perturbation is also conserved at third order on large scales in the adiabatic 
case. It is worth noting that higher order perturbation theory, as discussed in this 
chapter, has the advantage of being valid on all scales whereas fully non-linear methods, 
such as separate universe approaches are gradient expansions (in powers oik/ aH) , and 
so are only valid on superhorizon scales. 

Another application of our third order variables and equations, in particular the 
Klein-Gordon equation (5.49), is the calculation of the trispectrum by means of the field 
equations. Whereas calculations of the trispectrum so far derive the trispectrum from 
the fourth order action, it should also be possible to use the third order field equations 
instead. The equivalence of the two approaches for calculating the bispectrum, using 
the third order action or the second order field equations, has been shown in Ref. [136]. 
Having included tensor as well as scalar perturbations it will be in particular interesting 
to see and be an important consistency check for the theory whether we arrive at the 
same result as Ref. [137]. 

A final advantage of extending perturbation theory to third order is that, in doing 
so, one obtains a deeper insight into the second order theory. Also second order per- 
turbation theory, despite remaining challenging, becomes less daunting having explored 
some of the third order theory. 
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6.1. Summary 

The main focus of this thesis has been the study of cosmological perturbations beyond 
hnear order. In Section 1.1 we introduced the standard cosmological model, giving 
the background evolution and constraint equations and briefly discussing inflationary 
cosmology. In Chapter 2 we introduced the theory of cosmological perturbations up 
to second order, presenting the perturbed metric tensor and energy momentum tensor 
for both a perfect fluid, including all types of scalar, vector and tensor perturbations, 
and a scalar field. We then considered the behaviour of the perturbations under a 
gauge transformation in the active approach, using this behaviour to define gauges 
and construct gauge invariant variables. Next, we discussed the thermodynamics of a 
perfect fluid and defined the non-adiabatic pressure perturbation, closing the chapter by 
considering how non-adiabatic pressure perturbations can arise naturally in multiple- 
component systems. 

In Chapter 3 we continued the discussion of the foundations of cosmological per- 
turbation theory by presenting the dynamic and constraint equations, from energy 
momentum conservation and the Einstein field equations, up to second order in the 
perturbations. Starting with the linear theory, we gave the governing equations for 
scalar, vector and tensor perturbations of a perfect fluid in a gauge dependent form, i.e. 
without fixing a gauge. We then presented the gauge invariant form of the equations for 
three different gauges: the uniform density, the uniform curvature, and the longitudinal 
gauges, solving the equations for the latter two in the case of scalar perturbations. We 
then presented the Klein-Gordon equation for a scalar field, to linear order for both a 
canonical and non-canonical action, and highlighted the important difference between 
the adiabatic sound speed and the speed with which perturbations travel in a scalar 
field system (the phase speed). Finally, we investigated the perturbations of a system 
containing both dark energy and dark matter. 
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Having laid the foundations of the hnear order theory we then discussed the second 
order theory, presenting the governing equations for a perfect fluid derived, as at first 
order, from energy momentum conservation and the Einstein field equations. We then 
presented all equations in the uniform curvature gauge, which we then used in Chap- 
ter 4, and for scalars only in the Poisson gauge, in order to connect with the literature. 

In Chapter 4 we used the tools developed in the previous chapters to investigate 
non-linear vector perturbations in the early universe. Using the qualitative difference 
between the linear theory and the higher order theory we showed that, at second order in 
perturbation theory, vorticity is sourced by a coupling quadratic in linear energy density 
and entropy perturbations, extending Crocco's theorem from a classical framework to 
an expanding, cosmological background. In order to show this, we first defined the 
vorticity tensor in General Relativity and calculated the vorticity tensor at first and 
second order in cosmological perturbation theory using the metric tensor and fluid four 
velocity presented earlier. We then computed the evolution of the vorticity tensor by 
taking the time derivative, and using the governing equations from Chapter 3 to simplify 
the expressions and replace the metric perturbation variables. We found that at linear 
order the vorticity is not sourced in the absence of anisotropic stress, in agreement with 
the previous results known in the literature. However, at second order we obtained the 
novel result that there exists a non-zero source term for a fluid with a general equation 
of state (depending on both the energy density and entropy) which is quadratic in linear 
energy density and non-adiabatic pressure perturbations. 

Having derived this qualitative result we then gave the first quantitative solution, 
estimating the magnitude and scale dependence of the induced vorticity using simple 
input power spectra: the energy density derived in Chapter 3 approximated for small 
kr], and an ansatz for the non-adiabatic pressure perturbation. We found that the 
resulting spectrum has a surprisingly large magnitude, given that it is a second order 
effect, and a dependence on the wavenumber to the power of at least seven, given our 
assumptions. Thus, this spectrum is hugely amplified on small scales, rendering the 
vorticity not only possibly observable, but also important for the general understand- 
ing of the physical processes taking place in the early universe. 

In Chapter 5 we extended the formalism of cosmological perturbation theory from 
the second order theory to third order, starting with the gauge transformation rules 



6.2: Future Directions 



121 



and defining gauge invariant variables. Then, considering perfect fluids and scalar, 
vector and tensor perturbations we presented the energy and momentum conservation 
equations and the Klein-Gordon equation for a scalar field without fixing a gauge. 
Finally, we gave the components of the Einstein tensor at third order also in a gauge 
dependent form. 

6.2. Future Directions 

The work presented in this thesis can naturally be extended in several directions. One 
clear extension is to study the vorticity generation in specific models, moving beyond 
the simple ansatz for the non-adiabatic pressure perturbation used in Chapter 4. As 
discussed in Section 2.5.1, non-adiabatic pressure perturbations naturally arise in any 
system consisting of more than one component, such as a multiple fluid, or multiple 
scalar fields model. Even in the case of zero intrinsic non-adiabatic pressure perturba- 
tion, there exists a relative non-adiabatic pressure perturbation between the different 
components of the system which is proportional to the relative entropy perturbation, 
i.e. 



This can then be written in terms of field variables using the definitions in Chapter 2 
in order to obtain an expression for the relative entropy perturbation in multi-field 
inflation models. 

Another way in which this work can be extended is to exploit the potential for this 
mechanism to generate magnetic fields. As mentioned in Chapter 4, magnetic fields 
and vorticity are intimately related. This relationship has been studied in some detail 
in classical fluid mechanics and in astrophysical situations (see, e.g., Ref. [1 5.'^] for a 
review), though there is still much work to be done on incorporating magnetic fields into 
cosmological perturbation theory. Since the energy density is related to the magnetic 
field, 5*, through the expression 



when considering linear perturbations of the energy density one often considers pertur- 
bations of the magnetic field to 'half order'. That is, one assumes that the magnetic 
energy density as shown above to be of order 6^, is formally the same order as the 




(6.1) 
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scalar density perturbation, ensuring that the perturbed version of Eq. (6.1) holds at 
linear order (see, e.g., Ref. [26]). However, one does not need to use this technique and, 
in fact, when considering higher order perturbations it is not immediately clear how 
this will work. Instead, one can consider and develop cosmological perturbation theory 
to consistently include magnetic fields at integer order. Having done that, it will be 
possible to obtain estimates of the primordial magnetic field produced by the vorticity 
generated using cosmological perturbations as shown in this thesis by using the simple 
ansatz considered in Chapter 4 [ ] . 

A further extension to this work will be to consider the primordial magnetic field gen- 
erated from the relative non-adiabatic pressure perturbation in a system with multiple 
components, such as a hybrid inflation model. 

6.3. Outlook 

Cosmological perturbation theory has matured over the last few decades and has been 
incredibly successful in making predictions that agree with observations. However, with 
the data sets available to us continually growing in their size and quality it is now a 
realistic aim to use perturbation theory even beyond linear order to make predictions 
which are observationally testable. 

The main observable with which we can constrain our cosmological models is the 
CMB, of which we have, to date, collected much information, predominantly from 
the successful Wmap experiment. The Planck satellite [1] will greatly improve the 
temperature measurements of the CMB and together with the proposed CMBPoL 
satellite [21], will measure the polarisation of the background radiation. Since the 
CMB is not affected by the astrophysics of the late universe, usually one prefers the 
use of CMB data over other techniques for the study of higher order observables. 

However, with the recent technological advances. Large Scale Structure (LSS) sur- 
veys, such as the Sloan Digital Sky Survey (SDSS) [' 17] and the proposed 21cm 
anisotropy maps are attracting more attention as a way to probe the evolution of 
the universe at different epochs of its history. The 21cm signal, generated by neutral 
Hydrogen left over after the Big Bang, can probe the era after decoupling but before 
galaxy formation, i.e. between redshift 200 and 30, while LSS surveys probe out to 
around redshift 1. The 21cm anisotropy maps contain much more data than the CMB 
[95], though it should be noted that it is still not clear whether the foregrounds can be 
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removed with enough accuracy to enable reliable results. 

This is but one area where calculations at higher order can be tested against observa- 
tions. Thus, the future study of cosmological perturbation theory will greatly increase 
our understanding, serving to broaden and deepen our knowledge of the universe in 
which we live. 



A. Second Order i — j Einstein 
Equation 



In this appendix we give the second order {i — j) component of the Einstein Equation 
omitted from Section 3.2: 

Cijfe/ +C{ 2C5 j,fe+2'H(Cij7;,*+C5fej — CJj,fc)+| (i3ij/fe+-B]; jfc — 2i?ife/j) 

+{i?f-2Cf j+cJ,,'=+0i,''O(C'ijfe/+c;fc,j-q,,fc)+|i?5(Bife,/-v2Bi,+4W0i^ 



+2Cf 

-|(Bifc/Bf,j+-Bij,*=Bl,fc)+</'i 



(Bij/ +BI .+2<Ai/j+2W(Bij/+B5,j)-2C; ^-4WC; 



— ^IjZ, k—^il,jk-r^ij^kl 



+'5'j|(«^-^)(4</'?-Bifc^?f)+2(/,i 
+2Cf 



2C{',,fc-2C;fe;/+V2i?ife-Bi,,fc'+2W(B;^-<^i_fc-2C5fe,,+C{,,fc) 



+<^l'(Cl'''fc-Bl^fc)-|(2Bl^fcBl^/-Bl^,fcBl^'-3Bi,fcBl^,'^)+0l,fc(Cl^,'=-2Cl'^^-0l,'=)| 
=87rGa2|5P2<5"j+2(po+Po)i'iHfij+Bij)| • (A.l) 
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B. Third Order Einstein Tensor 



Here, we give the third order corrections to the components of the Einstein tensor. We 
do not spht up perturbations order by order. 



%=2C'-'' 2Cife(C'(j'=-2C'=\j7+V2C/-2«Cj*)+(2Cjfc,,-Cij,fe)(C^,*-2C^/)+Cii,fc(3C/,'=-C/,^ 



+2C^^B^ 



25* {C'l^ij—Cji^k) 



+Bi{2Cjk,''-C>'k,j)+Bicl>j+{2Bj^i-C'iM 



+mcj>'^{B\i-C-^\)+{Bi^j-Bj^i){\BiC''k,' -B'Ci^ ,k) 

Bi(2C,k,''-C''k,i+'t>,^-'2'^^B,)+B\^<i>-2C'\ 



+Bk,jBi{C''\i-2Ci'' ,J)+2HBi 



Ci'.Bi+2(t>(2Ciji-Cij,i+24>,i-^HBi) 



+BiBi{2Cjk,i''-C''kAj-^'^Cij)+B^Ci'.(C^k,i--2Cik!'~^,i^^^ (B.l) 



125 



B: Third Order Einstein Tensor 



126 



G°i=Ckj[2C^'^(c'',/-Cuy)+CiJ{c\^''-2C''\l-2<t>,'')+C'^i(C^''^^ 



+2{C'-Bf,j)<p,i+2<P{2C'^-2C'+Bi^,,-Bk,,j) 



+2Bi 



C''^C''i^ik—Cki,i'—Cu^k'+V'^Cik) 



+C'''{Cki,i^-C^i,ik-Cu,k'+Ci^M)MC'i,''-C''\i){Ch^i+Cik,^-Cii^k^^^^ 



+2q.</.(C"=fc/-2C"=J\fc)+B^ 



+Sj,i(B*,fe-C'='J+|Bj,fe(2Cf-S'=,iSi,'=) 



-20 



+C;^^C*J■,i+2C^■;.</.,i-2C^0/-Bfc,,■(Ci^'=-C'i'=,J■)-2B,•/>,i 



+2Cfc^C^'(Bi,,'=-B(,i'=) 



+B3 



^i(C'=fc,,-2Cjfc,*+4HB,)+2,/.(C*fc,„-C,fc,^-Cifc,/+V2Ci,+WB,-i)+^^ 



(Cjj^/j —Cj/j J -|-Cj/j_j) 



+2,^2(2Ci^._.-2q^./-Bj,iJ+V2Bi-12-K</.,i), 



(B.2) 
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+4Ci''Ci^{CkrnJl-Clrn,kj+CjrnM-CjkArn)+^Cl'^Ckl{Cj^^^^ 



m I J. Ij3ll 



B'0(Cfcij+C,^fc-C,fc,0+4C,>„,i(C'fc','"-C'fc'"/)-2Q^,,C'''-,fc 



-2Cj^(C';-B\,-0O+2C;,{2C';-Bfe/)+(:7';(Bfc,,+B,,fe)-2C^,B,-/-(</>'+B\O(^^ 
+S;jS',fc+S,-,B^/-2</.,fc0,,-2</.(s;,,.+B^.^+</.,,-^+4HCH(2C';.-S^j-S^^ 



2{C^rn,l—2ClrnJ"){C'^j,k — C'^k,3—Cjk,^)+2{2C'^l^k—CklJ^){Cjrn,^+C'^rn 



+'^{Cji,m—Cjm,i){C^^lP'—C^k,"^)—2C^'jC'f.i-\-'iC'''f.C'-i+{B'^^j+B^ 



+2C%Bi^k-2C%Bj^l-2C'.^B\i+Bk;Bij-2B'^{Cji;+CU,j-Cij,l) 



-Bf'Bk 



C%'-i(S*'.+Bj/') 



-4>;j-H{B\j+Bji-2Ci'.) 



+B*B'<=' 



C'i,jk-C^k,jl-C^j,kl+'^''Cjk+2HBk,. 



+B''B'- 



CklJ'j~CjlJ'k~C^l,jk+C'^j,kl 



C'ikiCh,*'-2C'^\l-'t>,*')+C'^lC''\i + \4',HBk,3+Bj,k) 



+(^^C'-i,>'-C>'\i){Bk,j-Bi^k)+Bi,k{C''j,'^-C'^j,^)-cl>,j{B\k-C'''^ 



+m<t>^{2C^'.-B\j-Bj;) 



+B 



C»fe{2C';_^-2C,,/'-B^^i+V2B,0+C"(4C;,_,-2q,_.-q,_,+2i?fe.,,-B;,,,-B,,fci) 

+C^k{4Ciri-2C\j-Cji;'+2Bi/-Bj/-B\ji)+Ci'j{2CMyB'^^^^ 

+2C'ji{C'\k-C'k,')+k(B'k+Ch,k-2Cu,%Bj;+B\j)+{C'j;+C'\^ 

+{Cjk;+C'kj-Cij,k){C'l-4>'-B\i)+{B\'-4HC'^){Cjk,i-^^^^ 



■\4>,k(B\3+Bj:)+Bj^i{c\y&\k)+Bk:4>,3+Bk,j<t>.' 



V2<t> 



B' C^-Cjk.^'+^C^'Bj-B 



\kj)^ 



+B>' \2Ci'.^^-C%.-Cjk,"+Bk/-\{Bj,k'+B\ik) j+iB'"+<P,''+4:nB>'){C'j^k-C'k,j-Cjk:) 
+C\{2C'''j-Bj^^)+2C^''C'^^-C'.^B\'' + \{B\j+Bj,%C^^^^ 
+40,j(2HB^-0/)-4«C<'=(2C';^.+Sfc,,+Bj.fc)+|(Bfc;S'=,,+B\fcBj,'=) 
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+C"=' 4C^„(C'"",H-2C(",fc'")+2q^(B'",(+Bi,'"-3C'";)+2(C"„,'"-2C'^",„)(2Q^,fc- 
H-C'^' +B'^'{4C^^fc-CM,rn)+2S^{2C,^™-C'"^,z)+C'"„,fc(4Q„"-C"„,0+C™n,fc{^^^ 

+CknM^Cl"r-'2Cr,n+'^Cln..k't','^-2C"^rn,k^,l-2Ckl,rn'^^^ 

+2(B,,fc-q,)(B-,^+</,'+4H0)+40(</,,H+B,',fc-q'()+Brri,fc{B(,'"-iS"',i) 



+2B. ( 2C,.,'=-2C'«;„,+B"^,.-V^B,-2H(b;+C'«.,-2C,.,'«-,,0) -iB,,.^ 

2C""'(C;,;-S;,^)+</.C';-8«02-4</.</.'-20B',,+B'(B('+C'"^,;-2C(^,'"-</.,(+2«Si) 



3fc ^C';'-B';(-V20-2H(2</.'+B',o) (| 



2Q^,fc'"-C""^,fc,-V2CH-2H(q,+B, 



C'"^-B''[^-V^<P+2H(H<p-3<P'-B''^k) 



+C'="0^2B;,fc-3C;;j+iSi,fc</.(S'=/-3B','=) 



+B\fc0(4</,'+B',O+20{Cfc,,'=<^/-C',,fc0.'=+20,fc</,,'=)+2B'='(2Cfc,/</.-0^^ 

+ B* C""m/(Bfe_;— B;_fe — 2C^j)+C''",m(4C^; — 2B(;j) + B'_;(2Cfcm,'" — C""m,fc — 

+Bm,KC'fc,'"-3C-fc/+C'-,fe)+q^(Cfc-/-C'-,fc)-</.'(C',,fc-2C'fei.')+B',feC^,- 
+20(^2Cfc(/'-2C";,,+Bi,fc'-V2Bfc-4?^(2B;;+C",,fc-2CH/-2</.,fc)) . 
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